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ABSTRACT 

Vascular  targeting  agents  (VTA)  are  new  types  of  anticancer  drugs  that  act  on  existing  tumor  vasculature, 
causing  vascular  disruption,  which  ultimately  leads  to  extensive  ischemic  tumor  cell  death.  One  major  goal  of 
this  projeet  is  to  assess  physiologieal  changes  at  different  time  points  following  VTA  in  breast  tumors. 
Continuous  studies  of  tumor  vascular  perfusion  and  tissue  oxygenation  by  MRI  confirmed  our  previous  findings 
during  Year  1.  Complementary  data  during  this  annum  showed  that  a  significant  drop  in  mean  tumor  p02  within 
90  min  after  administration  of  combretastatin  A4  phosphate  (CA4P)  and  a  further  decrease  was  observed  at  2  h. 
Intriguingly,  the  initial  changes  in  p02  in  central  and  peripheral  regions  were  parallel,  but  by  24  h  post  treatment 
significant  difference  was  apparent:  p02  in  the  periphery  improved  significantly,  while  the  center  remained 
hypoxic.  These  data  are  consistent  with  DCE  (dynamic  contrast  enhanced)  MRI,  which  revealed  a  -70% 
decrease  in  perfusion/permeability  at  2  h,  which  recovered  fully  after  24  h  in  a  thin  peripheral  region,  but  not  the 
tumor  center.  Based  on  the  imaging  results,  the  radiation  treatment  has  been  designed  and  initiated.  We  believe 
that  quantitative  p02  measurements  are  potentially  important  for  optimizing  therapeutic  combination  of  CA4P 
with  irradiation. _ 
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Background: 

Tumor  growth,  survival  and  metastasis  depend  critically  on  the  development  of  new  blood  vessels 
[1].  Therefore,  extensive  research  has  focused  on  developing  strategies  to  attack  tumor  vasculature  [1,2]. 
Tubulin  binding  agents,  e.g.,  combretastatin  A-4-phosphate  (CA4P)  represent  one  kind  of  vascular 
targeting  agent  (VTA)  [3,4].  Promising  preclinical  studies  have  shown  that  such  agents  selectively  cause 
tumor  vascular  shutdown  and  subsequently  trigger  a  cascade  of  tumor  cell  death  in  experimental  tumors 
[4,5].  Although  massive  necrosis  can  be  induced,  tumors  usually  regrow  from  a  thin  viable  rim.  Thus,  a 
combination  of  VTAs  with  additional  conventional  therapeutic  approaches  will  be  required  [6,7]. 
Indeed,  several  studies  involving  the  combination  of  the  VTAs  with  irradiation  or  chemotherapeutic 
agents  have  shown  enhanced  tumor  response  [8-10].  To  better  understand  the  mode  of  action,  and  hence, 
optimize  such  combinations,  in  vivo  imaging  approaches  have  been  initiated  to  monitor  physiological 
changes  resulting  from  VTA  administration  [11-13].  Dynamic  contrast  enhanced  (DCE)  MRI  based  on 
the  transport  properties  of  gadolinium-DTPA  (Gd-DTPA)  is  the  most  commonly  used  imaging  approach 
to  study  tumor  vascular  perfusion  and  permeability  [14,15].  For  combination  with  radiotherapy, 
measurement  of  tumor  oxygen  dynamics  will  be  especially  important  since  hypoxia  affects  radiation 
response.  By  applying  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for 
Dynamic  Oxygen  Mapping)  MRI  [16],  dynamic  tumor  oxygenation  can  be  monitored  following  the 
treatment  with  CA4P. 

Body: 

Continuing  from  Year  1,  tumor  physiological  changes  after  CA4P  injection  have  been  extensively 
studied  in  Year  2.  Consistent  with  previous  results  of  Year  1,  DCE  MRI  showed  a  ~70%  decrease, 
averaged  over  a  group  of  13762NF  rat  breast  tumors  (n  =  5),  in  perfusion/permeability  (initial  area  under 
signal-intensity  curve  (lAUC))  at  2  h.  The  lAUC  recovered  fully  after  24  h  in  a  thin  peripheral  region, 
but  not  the  tumor  center  (Fig.  1  and  Table  1)  [17]. 

Pre  CA4P2hr  CA4P24hr 


lAUC 
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Figure  1.  Dynamic  contrast  enhanced  (DCE)  MRI  with  respect  to  CA4P  (30  mg/kg)  treatment.  A) 
Normalized  T1 -weighted  contrast  enhanced  images  were  acquired  23  s  after  a  bolus  injection  of  contrast 
agent  pre,  and  2  h  and  24  h  after  treatment.  Significantly  less  signal  enhancement  was  observed  for  the 
whole  tumor  region  2  h  after  treatment.  A  full  recovery  was  apparent  in  the  tumor  rim  24  h  post¬ 
treatment.  B)  lAUC  frequency  of  DCE  MRI  on  voxel  by  voxel  basis  obtained  on  the  same  tumor  at  pre, 
2  h  and  24  h  after  CA4P.  Compared  with  pre,  curves  of  2  h  and  24  h  after  showed  a  significant  increase 
in  the  number  of  voxels  with  no  signal  enhancement  (typically  lAUC  <  0.05).  However,  the  frequency 
of  highly  enhancing  voxels  (lAUC  >  1 .5)  at  24  h  recovered  to  the  pre-treatment  level. 


Table  1  Normalized  lAUC  by  dynamic  contrast  enhanced  MRI 


Groups 

Mean  lAUC 

baseline 

2h 

24  h 

CA4P  (n  =  6) 

periphery 

0.88  ±  0.06 

0.31  ±0.08*^ 

0.86  ±  0.06 

(30  mg/kg) 

center 

0.44  ±  0.1  r 

0.13  ±0.08*^ 

0.12  ±0.05*'"^ 

Control  (n  =  3) 

periphery 

1.09  ±0.04 

0.93  ±0.13 

0.86  ±0.05 

(saline) 

center 

0.52  ±0.19^ 

0.46  ±  0.05^ 

0.44  ±  0.03^ 

LAUC;  initial  area  under  signal-intensity  curve; 

p  <  0.05  from  baseline;  ^p  <  0.05  from  periphery;  ^p  <  0.05  from  control. 


FREDOMMRl  further  demonstrated  that  a  significant  decrease  in  the  mean  pOa  for  the  group  of 
7  tumors  was  found  at  90  min  after  CA4P  (23  ±  5  vs.  9  ±  3  torr;  p  <  0.05)  and  a  further  decrease  at  2  h 
(2  ±  2  torr).  The  mean  pOa  increased  to  15  ±  4  torr  at  24  h,  which  was  still  significantly  lower  than 
pretreated  baseline  (p  <  0.05).  More  interestingly,  the  initial  changes  in  p02  in  central  and  peripheral 
regions  were  parallel,  but  by  24  h  post  treatment  significant  difference  was  apparent:  p02  in  the 
periphery  improved  significantly,  while  the  center  remained  hypoxic  (Figs.  2  and  3)  [17]. 


torr 
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Figure  2  A)  Distribution  of  hexafluorobenzene  (HFB)  in  a  representative  tumor.  An  overlay  of  '^F 
signal  density  on  the  anatomic  image  indicates  HFB  in  both  peripheral  and  central  regions.  B)  p02  maps 
obtained  from  the  same  tumor  comprise  two  separated  groups  of  voxels,  which  correspond  to  the 
locations  of  HFB  on  the  anatomic  image.  Twenty  five  individual  voxels  were  traceable  from  the 
pretreated  baseline  to  2  h  after  CA4P  with  oxygen  breathing.  Thirty  one  voxels  could  be  followed  24  h 
post-treatment.  Significant  decrease  in  p02was  evident  for  all  the  individual  voxels  after  CA4P,  and  p02 
did  not  respond  to  oxygen  inhalation  after  2  h.  The  24  h  maps  showed  improved  p02  and  significant 
response  to  oxygen  breathing  in  the  peripheral  region  (right),  but  not  in  the  central  region  (left). 


Figure  3  Mean  p02  curves  are  shown  for  the  peripheral  (■)  and  central  (o)  voxels  of  the  tumor  shown  in 
Fig.  2.  Significant  decrease  in  p02  was  found  as  early  as  30  min  after  CA4P  (30  mg/kg)  for  both 
peripheral  and  central  tumor.  *  p  <  0.05  from  baseline  air,  *  p  <0.05  from  24  h  air,  ^  p  <  0.05  from 
periphery. 

As  a  normal  tissue  control,  p02  in  femoral  muscle  was  also  monitored  before  and  after  injection  of 
CA4P.  Unlike  the  tumor  behavior,  normal  muscle  showed  no  decreased  p02  at  any  time  points  after 
CA4P  (Table  2)  [17,18].  These  results  suggest  that  the  vascular  targeting  agent,  CA4P  selectively  attack 
the  tumor  blood  vessels  with  little  effect  on  normal  vasculature. 


Table  2  Tumor  oxygen  dynamics  assessed  by  ”f  MRI  with  respect  to  CA4P  treatment 


Groups 

no. 

animals 

PO2  (mean  ±  se  torr) 

Time  course 

baseline  60  min  90  min  2  h  oxygen 

24  h  base  oxygen 

CA4P 
(30  mg/kg) 

Breast  tumor 

■ 

23  ±5  13  ±5  9  ±3*  2  ±2*  5  ±3* 

15±4  58±lU 

Saline 
(0.25  ml) 

3 

31  ±  10  38  ±  16  36  ±  16  35  ±14  79  ±23* 

39  ±15  108  ±26'' 

CA4P 

Normal  muscle  (30  mg/kg) 

3 

24  ±1  25  ±8  25  ±7  24  ±5  79  ±18* 

NA  NA 

*  p  <  0.05  from  baseline;  ^  p  <  0.05  from  24  h  air  (24  h  after  injection  of  CA4P);  NA:  no  measurement. 
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I  have  also  performed  extensive  studies  on  tumor  blood  vessels,  perfusion,  and  hypoxia  in  the  13762 
rat  breast  tumors.  Immunohistochemical  studies  of  tumor  hypoxia  and  vasculature  using  hypoxic  marker 
pimonidazole  and  endothelium  marker  CD31  showed  that  a  higher  microvascular  density  (MVD)  and 
lower  labeling  index  of  pimonidazole  (Figs  4  and  5). 


Figure  4  Comparison  of  pimonidazole  and  CD31  in  a  representative  tumor.  Fluorescent  immunostaing 
for  CD31  (left,  red)  and  hypoxia  (pimonidazole;  right,  green)  indicated  extensive  vasculature  and 
disperse  hypoxic  regions. 


CIXIl 


C.44P  2  h 


Cl-WP  24  h 


Figure  5  Perfusion  marker  Hoechst  staining  pre,  2  h,  and  24  h  post  CA4P  (left).  Vascular  endothelium 
of  the  same  field  was  immunostained  by  anti-CD31  (red,  middle).  A  good  match  (right)  between 
Hoechst  and  anti-CD31  stained  vascular  endothelium  was  found  in  the  pretreated  tumor.  Two  hours 
after  treatment,  significant  reduction  in  perfused  vessels  was  detected,  followed  by  a  recovery  at  24  h 
point.  Bar  =  100  pm. 
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As  proposed  time  frame,  the  Task  1  in  the  Statement  of  Work  is  complete.  One  peer-review 
paper  based  on  MRI  and  histological  studies  of  the  vascular  targeting  agent  has  been  published 
[17]. 

In  accordance  with  the  Statement  of  Work  Task  2,  experimental  radiation  therapy  has  been  designed 
and  initiated  based  on  the  MRI  oximetry  data.  Preliminary  data  of  control  group  (n  =  6)  and  irradiation 
alone  group  (n  =  6)  showed  a  single  dose  of  10  Gy  significantly  inhibited  tumor  growth  (Fig.  6).  Studies 
of  combination  of  radiation  and  CA4P  are  underway. 
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Days  after  irradiation 

Figure  6.  Tumor  growth  was  significantly  delayed  by  a  single  dose  of  10  Gy  radiation  (•)  compared  to 
control  tumors  (□). 


By  collaborating  with  Drs.  Hanli  Liu  and  Ralph  Mason  (Co-investigator),  non-invasive 
interrogation  of  breast  tumor  response  to  CA4P  has  also  been  performed  using  Near-Infrared 
Spectroscopy  (NIRS)  and  Bioluminescent  imaging  (BLI)  [19].  Consistent  and  complementary  results 
have  been  acquired  from  these  studies  (Figs.  7  and  8). 
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Figure  7  Qualitative  NIRS  detection  of  tumor  hemodynamic  changes  before,  2  h  (a)  and  24  h  (b)  after 
CA4P  (30  mg/kg;  i.p.)  administration.  Significant  decrease  in  concentration  of  both  oxygenated 
[Hb02]and  total  hemoglobin  ([Hb]totai;  blood  volume)  was  observed  immediately  after  CA4P.  In  contrast 
to  pretreated  baseline  study,  there  was  no  response  to  oxygen  inhalation  2  h  after  treatment.  Twenty  four 
hours  later,  significant  increase  in  [HbOj]  was  found  in  response  to  oxygen,  but  the  extent  of  increase 
was  about  6  times  lower  than  baseline  level. 
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Figure  8  Application  of  BLI  in  human  breast  tumor  xenograft.  A  nude  mouse  with  MDA-MB-231-luc 
mammary  tumor  was  monitored  temporally  before  and  after  treatment  with  the  vascular  targeting  agent 
combretastatin  (CA4P).  Significant  reduction  in  signal  intensity  was  observed  2  h  after  CA4P  (Right), 
compared  to  pre-treatment  (Left). 

Based  on  the  above  results,  one  peer  revievi^ed  paper,  one  paper  of  book  section  and  two 
abstracts  of  conference  proceeding  have  been  published  during  this  annum.  Several  other 
manuscripts  are  being  prepared. 

Key  Research  Accomplishments 

•  Assessment  of  dynamic  perfusion  and  oxygenation  in  the  breast  tumors  in  response  to  the 
vascular  targeting  agent,  Combretastatin  A4  phosphate,  by  in  vivo  MR  approaches 

a.  DCE  MRI  showed  significant  reduction  in  tumor  perfusion  and  permeability  2  h  after 
administration  of  CA4P  (30  mg/kg,  i.p.),  while  recovery  was  found  in  tumor  periphery  but  not 
center. 

b.  FREDOMMRl  revealed  that  tumor  p02  dropped  significantly  at  90  min  and  further  decreased 
at  2  h  after  CA4P. 

c.  Tumor  p02  improved  24  h  later,  but  it  was  still  significantly  lower  than  the  pretreated  baseline  p02 
With  respect  to  oxygen  breathing,  p02  in  the  periphery  improved  significantly,  while  the  center 
remained  hypoxic  pure  oxygen  inhalation  at  this  time  point. 

•  Correlation  of  MR  findings  with  biological  studies 

a.  Immunohistochemical  studies  of  tumor  hypoxia  and  vasculature  using  hypoxic  marker 
pimonidazole  and  endothelium  marker  CD3 1 . 

b.  Consistent  with  MRI  findings,  histological  study  of  tumor  perfusion  using  Hoechst  dye  33342 
showed  a  significant  reduction  in  perfused  vessels  at  2hr  after  CA4P,  which  recovered  24  h  later. 

•  Experimental  radiation  therapy 

a.  Irradiation  alone  on  tumor  growth  has  been  accomplished. 

•  Complementary  studies  using  other  non-invasive  techniques 

a.  Near-Infrared  Spectroscopy  (NIRS)  and  Bioluminescent  imaging  (BLI)  approaches  have  been 
applied  for  assessment  of  tumor  response  to  the  vascular  targeting  agent  CA4P.  The  data  is 
consistent  with  MRI  observation. 
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Reportable  Outcomes 

Reportable  outcomes  that  have  resulted  from  this  research  endeavor  in  this  annual  period  include; 

•  Peer  Reviewed  Publications: 

1.  Zhao,  D.,  Jiang,  L.,  Hahn,  E.W.,  and  Mason,  R.P.  Tumor  physiological  response  to  combretastatin 
A4  phosphate  assessed  by  MRI.  Int.  J.  Radiat.  Oncol.  Biol.  Phys.  62,  872-80, 2005. 

•  Sections  of  Edited  Books: 

1.  Kim,  J.  G.,  Zhao,  D.,  Mason,  R.  P.,  Liu  H.  Chemotherapeutic  (cyclophosphamide)  effects  on  rat 
breast  tumor  hemodynamics  monitored  by  multi-channel  NIRS.  SPIE,  5693, 282-92,  2005. 

•  Abstracts  (Published  Conference  Proceedings): 

1.  Zhao,  D.,  Jiang,  L.,  Adam,  A.,  Hahn,  E.W.,  and  Mason,  R.P.  Tumor  physiological  response  to 
antivascular  agent  combretastatin  A4  phosphate  assessed  by  magnetic  resonance  imaging.  96^* 
AACR,  #3799,  Anaheim,  CA,  Apr  2005. 

2.  Zhao,  D.,  Jiang,  L.,  Adam,  A.,  Hahn,  E.W.,  and  Mason,  R.P.  In  vivo  MRI  monitoring  of  breast 
tumor  response  to  the  vascular  targeting  agent  combretastatin  A4  phosphate.  Era  of  Hope  Meeting, 
P67-24,  Phaladephia,  PA,  June  2005. 

•  Manuscripts  in  preparation: 

Peer  reviewed  papers: 

1.  Zhao,  D.,  Jiang,  L.,  Hahn,  E.W.,  and  Mason,  R.P.  Evaluation  of  breast  tumor  microcirculation  and 

oxygenation  using  a  combination  of  BOLD,  DCE  and  '^F  MRI.  Magn.  Reson.  Med. 

2.  Zhao,  D.,  Harper,  A.,  Richer,  E.,  Liu,  Li.,  Slavine,  N.,  Shay,  J.  W.,  Antich,  P.  P.,  Mason,  R.  P.  Novel 

application  of  bioluminescent  imaging:  Interrogating  acute  effects  of  the  vascular  targeting  agent 
combretastatin. 

•  Employment  or  research  opportunity: 

Ya  Ren,  an  experienced  research  technician,  has  been  recently  recruited  as  a  Research  Assistant  by 
our  department.  Mrs.  Ren  will  contribute  40%  of  her  effort  to  this  project. 


Conclusion: 

Results  reported  here  were  successful  in  accordance  with  the  outlined  tasks  cited  in  the  original 
proposed  statement  of  work.  Five  publications  have  been  achieved  during  the  Year  2.  More  importantly, 
these  results  lay  foundation  for  the  undergoing  combination  treatment  of  vascular  targeting  agent  with 
irradiation  in  terms  of  order  and  timing.  I  am  confident  that  the  proposed  radiation  study  will  be  fulfilled 
by  the  next  term. 
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results  exhibit  a  linear  correlation  between  AHbO  and  ApOa  of  the  tumors 
under  hyperoxic  gas  intervention,  suggesting  that  the  NIRS  approach  could 
have  a  good  potential  value  in  the  clinic.  Finally,  the  newly  developed 
tumor  hemodynamic  model  allows  us  to  reveal  tumor  heterogeneities  at 
different  tumor  locations  based  on  the  multichannel  NIRS  results.  Through 
this  chapter,  we  lay  a  foundation  for  an  NIR  imaging  technique  to  be  fur¬ 
ther  developed  to  facilitate  investigations  of  tumor  heterogeneity  and  vas¬ 
cular  perfusion.  Such  a  noninvasive  imaging  approach  can  enhance  our 
understanding  of  the  dynamics  of  tumor  oxygenation  and  the  mechanism 
of  tumor  physiology  under  baseline  and  perturbed  conditions. 
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[18]  Measuring  Changes  in  Tumor  Oxygenation 

By  Dawen  Zhao,  Lan  Jiang,  and  Ralph  P.  Mason 

Introduction 

Significance  of  pO 2  in  Oncology 

It  has  long  been  appreciated  that  hypoxic  tumor  cells  are  more  resistant 
to  radiotherapy.'  Indeed,  a  3-fold  increase  in  radio  resistance  may  occur 
when  cells  are  irradiated  under  hypoxic  conditions  compared  with  oxygen 
pressure  p02  >  15  torr  for  a  given  single  radiation  dose.  However,  recent 
modeling  has  indicated  that  the  proportion  of  cells  in  the  range  0-20  torr 
may  be  most  significant  in  terms  of  surviving  a  course  of  fractionated  radio¬ 
therapy.^  Certain  chemotherapeutic  drugs  also  present  differential  efficacy, 
depending  on  hypoxia.^’'*  Increasingly,  there  is  evidence  that  hypoxia  also 


^  L.  Gray,  A.  Conger  et  al,  Br.  J.  Radiol.  26,  638  (1953). 

^  B.  G.  Wouters  and  J.  M.  Brown,  Radial.  Res.  147,  514  (1997). 
^  B.  Teicher,  J.  Lazo  et  al.  Cancer  Res.  41,  73  (1981). 

A.  C.  Sartorelli,  Cancer  Res.  48,  775  (1988). 
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influences  such  critical  characteristics  as  angiogenesis,  tumor  invasion,  and 
metastasis.^"^  Moreover,  repeated  bouts  of  intermittent  hypoxic  stress  may 
be  important  in  stimulating  tumor  progression.^  Thus  the  ability  to  mea¬ 
sure  p02  noninvasively  and  repeatedly,  with  respect  to  acute  or  chronic 
interventions,  becomes  increasingly  important. 

Early  work  examined  cells  in  vitro,  where  ambient  oxygen  concentra¬ 
tions  are  readily  controlled.  In  vivo,  hypoxia  may  be  achieved  by  clamping 
the  blood  supply  to  a  tumor,*°  but  other  levels  of  oxygenation  reflect  the 
interplay  of  supply  and  consumption.”’^^  Robust  fine-needle  polarographic 
electrodes  opened  the  possibility  of  measuring  p02  in  tumors  in  situ  and 
in  vivo  to  define  local  p02  under  baseline  conditions  or  with  respect  to 
interventions.  In  early  work.  Cater  and  Silver^^  showed  the  ability  to  moni¬ 
tor  PO2  at  individual  locations  in  patients’  tumors  with  respect  to  breathing 
oxygen.  Later,  Gatenby  et  al}'^  showed  that  p02  in  a  tumor  was  correlated 
with  clinical  outcome.  Tumor  oximetry  received  its  greatest  boost  with  the 
development  of  the  Eppendorf  Histograph  polarographic  needle  electrode 
system.^^  This  computer-controlled  device  equipped  with  a  stepper  motor 
can  reveal  distributions  of  tumor  oxygenation  and  has  been  applied  exten¬ 
sively  to  clinical  trials.  Many  reports  have  now  shown  that  tumors  are 
highly  heterogeneous  and  have  extensive  hypoxia;  furthermore,  strong 
correlations  have  been  shown  in  cervix  and  head  and  neck  tumors  between 
median  p02  or  hypoxic  fraction  and  survival  or  disease-free  survival.'^'^^"^® 
Extensive  hypoxia  also  has  been  found  in  tumors  of  the  prostate  and 
breast. Thus  tumor  oxygenation  is  now  recognized  as  a  strong 


^  E.  K.  Rofstad,  K.  Sundfor  et  al.,  Br.  J.  Cancer  83,  354  (2000). 

®  K.  De  Jaeger,  M.  C.  Kavanagh  et  at.,  Br.  J.  Cancer  84,  1280  (2001). 

’  M.  Hockel  and  P.  Vaupel,  J.  Natl.  Cancer  Inst.  93,  266  (2001). 

®  H.  J.  Knowles  and  A.  L.  Harris,  Breast  Cancer  Res.  3,  318  (2001). 

®  R.  A.  Cairns,  T.  Kalliomaki  et  al.  Cancer  Res.  61,  8903  (2001). 

J.  Moulder  and  S.  Rockwell,  Int.  J.  Radial  Oncol  Biol  Phys.  10,  695  (1984). 

”  T.  W.  Secomb,  R.  Hsu  et  al,  Adv.  Exp.  Med.  Biol  454,  629  (1998). 

M.  W.  Dewhirst,  B.  Klitzman  et  al,  Int.  J.  Cancer  90, 237  (2000). 

D.  Cater  and  I.  Silver,  Acta  Radiol  53,  233  (1960). 

^''R.  a.  Gatenby,  H.  B.  Kessler  et  al.  Ini  J.  Radiat.  Oncol  Biol  Phys.  14,  831  (1988). 
M.  Nozue,  I.  Lee  et  al,  J.  Surg.  Oncol.  66,  30  (1997). 

D.  M.  Brizel,  S.  P.  Scully  et  al.  Cancer  Res.  56,  941  (1996). 

M.  Hockel,  K.  Schlenger  et  al.  Cancer  Res.  56,  4509  (1996). 

M.  Nordsmark,  M.  Overgaard  et  al,  Radiother.  Oncol  41,  31  (1996). 

A.  W.  Fyles,  M.  Milosevic  et  al,  Radiother.  Oncol.  48, 149  (1998). 

H.  Knocke,  H.  D.  Weitmann  et  al,  Radiother.  Oncol  53,  99  (1999). 
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P.  Hohenberger,  C.  Felger  et  al.  Breast  Cancer  Res.  Treal  48,  97  (1998). 

B.  Movsas,  J.  D.  Chapman  et  al.  Urology  53, 11  (1999). 


380 


GENERAL  METHODS 


[18] 


prognostic  indicator,  and  this  device  has  laid  a  convincing  foundation 
for  the  value  of  measuring  p02  in  patients.  However,  the  Histograph  is 
highly  invasive,  and  it  is  not  possible  to  make  repeated  measurements  at  in¬ 
dividual  locations,  precluding  dynamic  studies  to  assess  the  influence  of 
interventions  on  tumor  p02. 

Given  that  hypoxic  tumors  are  more  resistant  to  certain  therapies,  it  be¬ 
comes  important  to  assess  tumor  oxygenation  as  part  of  therapeutic  plan¬ 
ning.  Patients  could  be  stratified  according  to  baseline  hypoxia  to  receive 
adjuvant  interventions  designed  to  modulate  p02,  or  more  intense  therapy 
as  facilitated  by  intensity  modulated  radiation  therapy  (IMRT).  Tumors, 
which  do  not  respond  to  interventions,  may  be  ideal  candidates  for  hyp¬ 
oxia-selective  cytotoxins  (e.g.,  tirapazamine^'*).  Noting  that  any  therapy 
and  intervention  may  have  side  effects  or  simply  add  to  clinical  costs,  it  is 
vital  that  efficacy  be  established  and  therapy  be  optimized  for  an  individual 
patient.  Whether  initially  hypoxic  regions  of  a  tumor  can  be  modified  to 
become  better  oxygenated  has  long  been  considered  a  key  to  improving 
outcome  of  irradiation.  However,  many  attempts  to  improve  therapeutic 
outcome  by  manipulation  of  tumor  oxygenation  have  shown  only  modest 
success  in  the  clinic,^^  and  it  is  thought  that  lack  of  success  may  have 
resulted  from  inability  to  identify  those  patients  who  would  benefit  from 
adjuvant  interventions. 

Although  p02  determinations  could  be  of  great  clinical  value,  they  are 
also  vital  to  many  laboratory  investigations  of  new  drugs  and  studies  of 
tumor  development.  Given  the  potential  importance  of  measuring  pOa, 
many  diverse  techniques  have  been  developed,  as  reviewed  by  others 
previously and  here,  in  the  next  section. 

Methods  of  Measuring  Tumor  Oxygenation 

Table  I  lists  various  techniques  that  have  been  reported  to  provide 
quantitative  estimates  of  p02.  Historically,  polarographic  needle  oxygen 
electrodes  have  been  considered  a  “gold  standard,”  and  they  have  been  ap¬ 
plied  in  the  clinic  since  the  1950s.  One  or  more  electrodes  may  be  placed  in 
a  tumor,  facilitating  measurement  of  baseline  p02  and  dynamic  response  to 


J.  M.  Brown,  Mol.  Med.  Today  6,  157  (2000). 

J.  Overgaard  and  M.  R.  Horsman,  Semin.  Radial.  Oncol.  6,  10  (1996). 

“H.  B.  Stone,  J.  M.  Brown  et  at.  Radial.  Res.  136,  422  (1993). 

R.  P.  Mason,  S.  Ran  el  al,  J.  Cell  Biochem.  87S,  (2002). 

H.  M.  Swartz,  Biochem.  Soc.  Trans.  30,  248  (2002). 

M.  Swartz  and  J.  F.  Dunn,  in  “Oxygen  Transport  to  Tissue  XXIV”  (J.  F.  Dunn  and 
H.  M.  Swartz,  eds.),  Vol.  530,  p.  1.  Kluwer  Academic,  New  York,  2003. 
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interventions.^^’^'’"'*^  Initially,  the  focus  was  on  generating  finer  needles, 
which  would  be  less  invasive,  and  tips  as  fine  as  a  few  microns  have  been 
applied  to  animal  tissues.^^  However,  such  needles  are  progressively  brittle 
and  generate  such  small  current  that  stray  electromagnetic  fields  can  inter¬ 
fere.  Stationary  electrodes  sample  limited  volumes,  and  recognizing  tumor 
heterogeneity,  the  Eppendorf  Histograph  was  developed  to  generate  mul¬ 
tiple  measurements  along  tracks  in  tumors.^^’^"*’^^  Following  extensive  stud¬ 
ies  in  animals,  the  Histograph  has  found  widespread  application  in  the 
clinical  setting  and  has  unequivocally  revealed  hypoxia  in  many  tumor 
types,  for  example,  head  and  neck,^®’"*’  cervix,^®’^®  breast, and  pros- 
tate.^^  Moreover,  p02  distributions  have  been  found  to  have  prognostic 
value.  Disease-free  survival  is  significantly  worse  for  patients  with  hypoxic 
tumors,  though  the  optimal  prognostic  parameter  has  variously  been 
median  p02  or  percent  measurements  <5  torr  (HF5). 

Although  the  Eppendorf  Histograph  uses  a  large  invasive  needle  (size  = 
26  G  or  about  0.35  mm),  it  has  provided  great  impetus  for  further  investiga¬ 
tions.  One  aspect  is  the  application  of  less  invasive  probes.  Fiber-optic 
probes  are  typically  finer  and  do  not  consume  oxygen  during  measurement. 
Typically,  only  two  or  four  locations  are  sampled  simultaneously,  but  as  with 
the  earlier  electrodes,  these  optical  probes  facilitate  observation  of  dynamic 
changes  in  p02  in  response  to  interventions.^^"^^  Both  the  current  commer¬ 
cial  systems,  the  OxyLite  (http://www.oxford-optronix.com/tissmon/oxylite/ 
oxylite.htm)  and  FOXY  (http://www.oceanoptics.com/Products/foxyfaqs. 
asp),  exploit  the  fluorescent  quenching  by  oxygen  of  a  ruthenium  complex 
coating.  OxyLite  measures  fluorescent  lifetime,  whereas  FOXY  uses  a 
simple  intensity  integration  and  is  correspondingly  much  cheaper.  Fibers 
are  fragile,  and  coatings  have  a  limited  lifetime. 
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Reporter  molecules  have  been  developed  for  use  with  electron  spin 
resonance  (ESR  or  EPR),  where  the  line  width  is  highly  sensitive  to 
oxygen?*’'*"^"'^^  Two  primary  approaches  are  used:  (1)  direct  intratumoral 
(IT)  injection  of  char  crystals, phthalocyanine,^’  or  India  ink’^**  into  a 
tissue  or  (2)  intravenous  (IV)  infusion  of  water-soluble  agents  which  dis¬ 
seminate  throughout  the  tumor  vasculature.'*^’^^  Direct  IT  injection  is  inva¬ 
sive  and  has  generally  been  applied  as  a  spectroscopic  approach  to  report 
pOa  at  single  locations  only.  Nonetheless,  significant  data  have  been 
achieved  demonstrating  hypoxiation  and  reoxygenation  with  respect  to  ir¬ 
radiation,  and  the  importance  of  timing  successive  radiation  doses  to  coin¬ 
cide  with  reoxygenation.^^  Char  particles  may  be  stable  in  tissue  for  weeks 
to  years,  allowing  measurements  of  chronic  changes  in  tissues  (e.g.,  accom¬ 
panying  tumor  growth).^*^  The  IV  approach  is  noninvasive,  but  reporter 
molecules  may  predominately  distribute  in  the  well-perfused  vasculature, 
potentially  biasing  measurements  toward  the  well-oxygenated  tumor 
regions.  Progressive  uptake  and  clearance  of  agents  produces  variable  con¬ 
centrations,  and  some  agents  degrade  in  tissue  requiring  appropriate  cor¬ 
rection  factors.'^^  Nonetheless,  images  of  tumor  oxygen  distribution  have 
been  reported,  including  three-dimensional  representations.'*’^  Spin  radicals 
also  may  be  applied  to  a  combined  ESR-NMR  (nuclear  magnetic  res¬ 
onance)  approach,  Overhauser-enhanced  magnetic  resonance  imaging 
(OMRI),  exploiting  the  Overhauser  enhancement  in  the  tissue  water 
proton  MRI  signal  that  occurs  by  polarization  transfer  from  free  radicals 
upon  electromagnetic  irradiation.^® 
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Vascular  oxygenation  has  been  probed  by  fluorescence  or  phosphores¬ 
cence  imaging  based  on  reporter  complexes  delivered  Historically, 

the  approach  was  limited  to  superficial  tissues  due  to  limited  light  penetra¬ 
tion.  The  latest  molecules  are  active  in  the  near-infrared,  permitting 
greater  depth  of  signal  penetration.® 

NMR  facilitates  interrogation  of  deep  tissues  noninvasively,  and 
NMR  approaches  will  be  reviewed  in  detail  in  the  following  section.  The 
methods  discussed  earlier  provide  direct  quantitative  measurements  of 
p02  based  on  various  physiochemical  parameters,  such  as  electric  current, 
fluorescent  lifetime,  magnetic  resonance  linewidth,  or  relaxation.  Other 
approaches  are  less  direct,  but  can  reveal  hypoxia  or  correlates  of  p02. 

Specific  classes  of  reporter  molecules  have  been  developed  to  reveal 
hypoxia^^’®  (e.g.,  pimonidazole,®’®  EF5,®’®  CC1-103F,®  Cu-ATSM'^°’’^ 
galactopyranoside  IAZA’^).  Following  IV  infusion,  these  agents  become 
reduced  in  tissues  and  are  trapped.  However,  in  the  presence  of  oxygen 
they  are  reoxidized  and  ultimately  clear  from  the  body.  Histologic  assess¬ 
ment  of  the  distribution  of  these  agents  provides  microscopic  indications 
of  local  hypoxia.  EF5,  pimonidazole,  and  Cu-ATSM  are  currently  being 
tested  in  clinical  trials,  and  correlations  have  been  reported  with  clinical 
'  outcome.^®’®’’®  Many  variants  have  been  proposed  over  the  past  20  years, 
and  incorporation  of  radionuclides  has  facilitated  noninvasive  investiga¬ 
tions  using  positron  emission  tomography  (PET)  or  single  photon  emission 
computed  tomography  (SPECT),  while  ^^F  labels  permitted  NMR  spec¬ 
troscopy.’^-’'^  Generally,  only  a  single  time  point  is  investigated,  but  dy¬ 
namic  variations  in  hypoxia  may  be  assessed,  even  in  biopsy  specimens, 
by  applying  pairs  of  hypoxia  reporters  in  a  pulse-chase  fashion  with  respect 
to  an  intervention,  as  shown  by  Ljungkvist  et 
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Several  studies  have  shown  a  lack  of  correlation  between  hypoxic 
marker  binding  and  pOa  assessed  using  the  Eppendorf  Histograph,  which 
may  be  related  to  chronic  versus  acute  hypoxia,  or  the  extent  of  necro- 
sis7^”’^  The  ultimate  value  of  the  techniques  is  evidenced  by  correlations 
between  uptake  and  outcome.^’™’’^’’**  Recent  data  also  indicate  that  EF5 
fluorescence  may  be  correlated  with  p02.^^ 

The  techniques  discussed  so  far  all  depend  on  exogenous  reporter  mol¬ 
ecules  or  probes.  Ideally,  oxygenation  could  be  related  to  endogenous 
characteristics.  Because  many  biochemical  pathways  are  under  oxygen 
regulation,  they  can  provide  an  elegant  window  on  hypoxia,  for  example, 
induction  of  hypoxia-inducible  factor  1  (HIF-1)  and  glucose  transporter  1 
(Glut-1)  together  with  secondary  responses,  such  as  increased  production 
of  vascular  endothelial  growth  factor  (VEGF),  N1P3  and  tumor-associated 
macrophage  activity.*  Such  molecules  indicate  hypoxia,  though  they  may 
be  induced  by  other  factors.  Intrinsic  radiation  sensitivity  also  may  be 
assessed  using  the  Comet  assay These  assays  each  require  biopsy.  Other 
markers  potentially  associated  with  hypoxia  may  be  found  in  the  plasma  or 
urine  and  have  been  correlated  with  clinical  outcome.*^  An  attractive  alter¬ 
native  is  the  introduction  of  transgenes  with  hypoxic  response  elements 
(HREs)  as  promoter  sequences  coupled  to  reporter  genes  such  as  GFP 
(green  fluorescent  protein)*^’**  or  luciferase.*'*’*”’  GFP  synthesis  is  an  ener¬ 
getic  process,  which  could  be  hindered  under  hypoxia  conditions.  Likewise, 
bioluminescence  accompanying  action  of  luciferase  on  luciferin  requires 
adenosine  triphosphate  (ATP)  and  O2,  but  reports  suggest  that  even  under 
exceedingly  low  p02,  sufficient  oxygen  remains  to  reveal  hypoxia. 

Many  practical  considerations  govern  clinical  application  of  oximetry 
methods.  Proton  MRI  is  routinely  applied  for  anatomic  evaluation  of 
tumors  and  would  provide  an  ideal  conduit  for  prognostic  investigations. 
Application  of  contrast  agents  may  reveal  tumor  boundaries  to  enhance 
detectability,  and  the  dynamic  contrast  enhancement  (DCE)  changes 
provide  insight  into  vascular  perfusion  and  surface  permeability  area.*® 
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Specific  studies  have  shown  a  correlation  between  DCE  and  p02,**’  and 
indeed,  a  theoretical  underpinning  has  been  provided  based  on  the  Krogh 
cylinder  model.®*^  However,  the  correlation  is  unlikely  to  be  widely  applic¬ 
able,  since  DCE  is  sensitive  to  vascular  flow,  perfusion,  and  permeability, 
where  pOa  depends  on  oxygen  consumption  as  well  as  delivery. 

Blood  oxygen  level  dependent  (BOLD)  contrast  proton  NMR  facili¬ 
tates  rapid  interrogation  of  vascular  oxygenation  and  is  particularly 
appropriate  for  examining  dynamic  responses  to  interventions.'*®’®®”®^ 
Deoxyhemoglobin  is  paramagnetic  and  induces  signal  loss  in  T2*-weighted 
images.  However,  BOLD  does  not  provide  absolute  p02  values  and  is  con¬ 
founded  by  the  influence  of  blood  flow,  as  investigated  extensively  by 
Howe  et  who  termed  the  expression  FLOOD  (flow  and  oxygen  level 
dependent)  contrast.  In  addition,  variation  in  vascular  volume  can  intro¬ 
duce  signal  perturbation.®^  Nonetheless,  some  studies  have  indicated  a 
correlation  with  relative  p02,  but  poor  indication  of  absolute  p02.^^’®'*’®^ 

Near-infrared  spectroscopy  (NIRS)  offers  an  alternative  approach 
based  on  the  differential  light  absorption  of  the  strong  chromophores  oxy¬ 
hemoglobin  and  deoxyhemoglobin.  NIRS  provides  a  noninvasive  means  to 
monitor  global  tumor  vascular  oxygenation  in  real  time  based  on  endoge¬ 
nous  molecules.  Although  many  NIRS  investigations  have  been  conducted 
in  the  brain  and  breast  in  both  laboratory  and  clinical  settings  over  the 
past  decade,  there  have  been  relatively  few  reports  regarding  solid 
tumors.'*^’®^"*”^  Most  studies  to  date  have  used  reflectance  mode.  By 
contrast,  we  have  favored  transmission  mode,  so  as  to  interrogate  deep 
tumor  regions,  and  we  have  presented  preliminary  studies  in  rat  breast 

R.  A.  Cooper,  B.  M.  Carrington  et  al,  Radiother.  Oncol.  57,  53  (2000). 

Z.  Wang,  M.-Y.  Su  et  al,  TechnoL  Cancer  Res.  Treat.  1,  29  (2002). 

A.  Al-Hallaq,  J.  N.  River  et  al,  Int.  J.  Radiat.  Oncol  Biol  Phys.  41,  151  (1998). 

S.  P.  Robinson,  F.  A.  Howe  et  al,  Semin.  Radiat.  Oncol  8,  198  (1998). 

R.  Mazurchuk,  R.  Zhou  et  al,  Magn.  Reson.  Imaging  17,  537  (1999). 

®^F.  A.  Howe,  S.  P.  Robinson  et  al,  Magn.  Reson.  Imaging  17,  1307  (1999). 

^  F.  A.  Howe,  S.  P.  Robinson  et  al,  NMR  Biomed.  14,  497  (2001). 

C.  Baudelet  and  G.  Gallez,  Magn.  Reson.  Med.  48,  980  (2002). 

X.  Fan,  J.  N.  River  et  al,  Int.  J.  Radiat.  Oncol  Biol  Phys.  54,  1202  (2002). 

®*H.  D.  Sostman,  S.  Rockwell  et  al,  Magn.  Reson.  Med.  20,  253  (1991). 

^  R.  G.  Steen,  K.  Kitagishi  et  al,  J.  Neurooncol  22,  209  (1994). 

®®F.  Steinberg,  H.  J.  Rohrbom  et  al,  Adv.  Exp.  Med.  Biol  428,  553  (1997). 

L.  Hull,  D.  L.  Conover  et  al,  Br.  J.  Cancer  79, 1709  (1998). 

B.  J.  P.  van  der  Sanden,  A.  Heerschap  et  al,  Magn.  Reson.  Med.  42,  490  (1999). 

H.  Liu,  Y.  Song  et  al,  Appl  Op.  39,  5231  (2000). 

M.  Kragh,  B.  Quistorff  et  al,  Eur.  J.  Cancer  37,  924  (2001). 

G.  Gulsen,  H.  Yu  et  al,  Technol  Cancer  Res.  Treat.  1,  497  (2002). 

E.  L.  Heffer  and  S.  Fantini,  Appl.  Opt.  41,  3827  (2002). 

T.  O.  McBride,  B.  W.  Pogue  et  al,  J.  Biomed.  Opt.  7,  72  (2002). 

J.  G.  Kim,  D.  Zhao  et  al,  J.  Biomed.  Optics  8,  53  (2003). 


390 


GENERAL  METHODS 


[18] 


and  prostate  tumors  with  respect  to  various  interventions."'^’^®^'^'’®  NIR 
approaches  are  presented  in  detail  in  Chapter  17  of  this  volume. 

Each  technique  has  specific  virtues  and  drawbacks,  which  must  be  con¬ 
sidered  for  any  given  application,  particularly  the  degree  of  invasiveness, 
the  ability  to  generate  maps  of  heterogeneity,  and  the  ability  to  assess  dy¬ 
namic  changes.  In  addition,  the  location  of  a  measurement  (e.g.,  vascular 
versus  tissue  compartments,  the  precision  of  measurements,  and  spatial 
and  temporal  resolution)  must  be  considered.  For  further  details  of  the 
techniques  described  earlier,  the  reader  is  referred  to  the  references.  In 
the  next  section,  we  present  NMR  approaches  in  greater  detail. 

NMR  Approaches  to  Measuring  p02 

Nuclear  magnetic  resonance  (NMR)  is  attractive  because  it  is  inherently 
noninvasive.  Liquid-state  NMR  is  characterized  by  several  parameters,  in¬ 
cluding  signal  amplitude,  chemical  shift  {6),  spin-spin  relaxation  (T2),  and 
spin-lattice  relaxation  (Ti).  Oxygen  could  be  quantified  using  NMR, 
but  this  is  rather  esoteric.'”’  Alternatively,  it  has  long  been  recognized  that 
the  oxygen  molecule  (O2)  is  paramagnetic,  causing  increased  spin-lattice  re¬ 
laxation  rates  (Ri  =  1/Ti).  Indeed,  physical  and  theoretical  chemists  must  go 
to  great  lengths  to  rigorously  remove  oxygen  from  solutions  (using  freeze- 
thaw  procedures)  to  achieve  inherent  relaxation  rates  for  studying  nuclear 
interactions.'”''  Proton  NMR  studies  have  reported  changes  in  the  water  re¬ 
laxation  rate  as  a  result  of  tissue  oxygenation,'””  but  many  other  processes 
(metal  ions,  cellularity,  pH,  ionic  strength)  also  cause  relaxation,  and  thus 
it  is  not  suitable  for  detecting  p02,  except  under  rare  circumstances,  such 
as  with  the  eye."”  There  is  also  a  substantial  temperature  response,  whereas 
the  relaxivity  due  to  oxygen  is  only  0.0002  s“'/torr."” 

However,  several  investigators  showed  that  the  '”F  NMR  spin-lattice 
relaxation  rates  for  fluorocarbons  are  much  more  sensitive  to  p02.'"’"’ 
Thomas  et  <2/."^  pioneered  the  application  of  '”F  NMR  relaxometry  to 
measure  p02  in  tissues,  in  vivo,  including  lung,  liver,  and  spleen;  several 
other  investigators  demonstrated  feasibility  and  applications,"'*"'’”  as 
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reviewed  some  years  ago  by  Mason.^^^  The  NMR  Ri  of  perfluorocar- 
bons  (PFCs)  varies  linearly  with  and  each  resonance  is  sensitive 

to  p02,  temperature,  and  magnetic  field,  but  importantly,  is  essentially  un¬ 
responsive  to  pFl,  CO2,  charged  paramagnetic  ions,  mixing  with  blood,  or 
emulsification.’^^"^^’ 

A  particular  PFC  molecule  may  have  multiple  resonances,  and  each 
resonance  has  a  characteristic  Ri  response  to  p02.  This  is  attributed  to 
steric  effects  of  O2,  as  it  approaches  the  molecule,^^^  which  implies  that 
perfluorinated  groups,  which  are  both  geometrically  and  magnetically  com¬ 
parable,  should  have  similar  Ri  responses  to  oxygen  tension.  At  a  fixed 
temperature  and  magnetic-field  strength,  the  Ri  response  to  p02  of  any 
single  resonance  obeys  the  simple  formula 

Ri  -  Ria  -F  (RipX)  (1) 

where  X  is  the  mole  fraction  of  O2  dissolved  in  the  PFC,  Rja  is  the  anoxic 
relaxation  rate,  and  Rip  is  the  relaxation  rate  due  to  the  paramagnetic  con¬ 
tribution  of  oxygen.  According  to  Henry’s  law,  the  dissolved  mole  fraction 
is  related  directly  to  the  partial  pressure  of  oxygen, 

PO2  =  KX  (2) 

where  K  represents  Henry’s  constant  for  a  given  solution  of  gas  at  a  speci¬ 
fied  temperature.  By  substitution, 

Ri  =  Ria  +  (Rip/K)p02  (3) 

The  slope  (RiplK)  indicates  the  response  of  a  particular  resonance  to  p02. 

PFCs  essentially  act  as  molecular  amplifiers,  since  the  solubility  of 
oxygen  is  greater  than  in  water,  but  thermodynamics  require  that  the  p02 
in  the  PFC  will  rapidly  equilibrate  with  the  surrounding  medium,  and 
estimates  of  diffusion  suggest  the  equilibration  can  occur  within  seconds. 
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Because  relaxation  is  proportional  to  oxygen  concentration,  the  effect 
will  be  greater  at  a  given  p02  than  for  water.  Importantly,  ions  do  not  enter 
the  hydrophobic  PFC  phase,  and  thus  do  not  affect  the  bulk  relaxation. 
Indeed,  PFCs  are  typically  exceedingly  hydrophobic  and  do  not  mix 
with  the  aqueous  phases,  but  rather  form  droplets  or  emulsions.  Based  on 
these  principles,  PFCs  have  been  applied  to  in  vivo  p02  measurements. 
Characteristics  of  many  diverse  PFCs  are  summarized  in  Table  II. 

At  any  given  magnetic  field  (Bo)  and  temperature  (T),  sensitivity  to 
changes  in  p02  is  given  by  Ri  =  a  +  bp02.  Thus  a  greater  slope  is  impor¬ 
tant,  and  the  ratio  77  =  b/a  has  been  proposed  as  a  sensitivity  index.^^*^  Gen¬ 
erally,  a  small  “a”  value  (intercept)  represents  greater  sensitivity,  but  it 
also  generates  longer  Ti  values  under  hypoxic  conditions,  potentially  in¬ 
creasing  data  acquisition  times.  Indeed,  the  Tj  of  hexafluorobenzene 
(HFB)  at  4.7  T  may  reach  12  s,  potentially  creating  long  imaging  cycles, 
but  this  is  readily  overcome  by  applying  single-shot  (echo  planar)  imaging 
techniques,  as  presented  in  a  later  section. 

Many  PFCs,  such  as  perfluorotributylamine  (PFTB),  perflubron  (for¬ 
merly  referred  to  as  perfluorooctyl  bromide;  PFOB),  and  Therox  (F44- 
E),  have  several  NMR  resonances,  which  can  be  exploited  to  provide 
additional  information  in  spectroscopic  studies,  but  seriously  hamper  ef¬ 
fective  imaging.  Multiple  resonances  can  lead  to  chemical  shift  artifacts 
in  images,  which  compromise  the  integrity  of  relaxation  time  measure¬ 
ments,  though  they  can  be  avoided  by  selective  excitation,  or  detection, 
chemical  shift  imaging,  deconvolution,  or  sophisticated  tricks  of  NMR 
spin  physics."®’^’®’'^®"^^^  These  approaches  add  to  experimental  com¬ 
plexity  and  are  generally  associated  with  lost  signal  to  noise  ratio  (SNR). 
Thus  we  strongly  favor  PFCs  with  a  single  resonance,  and  we  will  describe 
the  use  of  though  some  research  groups  favor 

5-crown-5-ether 


S.  R.  Thomas,  R.  G.  Pratt  et  al,  Magn.  Reson.  Imaging  14,  103  (1996). 

L.  J.  Basse,  R.  G,  Pratt  et  al.,J.  Comp.  Ast.  Tomagr.  12,  824  (1988). 

R.  P.  Mason,  N.  Bansal  et  al,  Magn.  Reson.  Imaging  8,  729  (1990). 

H.  K.  Lee  and  O.  Nalcioglu,  J.  Magn.  Reson.  Imaging  2,  53  (1992). 

U.  Noth,  R.  Deichmann  et  al,J.  Magn.  Reson.  B  105, 233  (1994). 

R.  G.  Pratt,  J.  Zheng  et  al,  Magn.  Reson.  Med.  37,  307  (1997). 

D.  Zhao,  A.  Constantinescu  et  al,  Radiat.  Res.  156,  510  (2001). 

Y.  Song,  A.  Constantinescu  et  al,  Technol  Cancer  Res.  Treal  1,  471  (2002). 

D.  Zhao,  A.  Constantinescu  et  al,  Radiat  Res.  159,  621  (2003). 

D.  Zhao,  S.  Ran  et  al.  Neoplasia  5,  308  (2003). 

^^®B.  P.  J.  van  der  Sanden,  A.  Heerschap  et  al,  Int  J.  Radiat  Oncol  Biol  Phys.  44,  649 
(1999). 

T.  Q.  Duong,  C.  Ladecola  et  al,  Magn.  Reson.  Med.  45,  61  (2001). 


TABLE  II 

NMR  Characteristics  and  Applications  of  PFCs  for  Tissue  Oximetry 


[18] 


DYNAMIC  ASSESSMENT  OF  TUMOR  OXYGENATION 


393 


r.--  3  ® 

f  e  ^ 

c  S  <u  g 
«  g  c  S 
»i  o  J  4  ' 
c  ^  ^ 

S  M  o  C 
«  S  Jc 
S  eg  § 
K 


Vn  TJ  -T 

■S 

C  >  c 


.  .  B 

tf  tH  O 

1  i  &• 

si  l« 

“  s  s? 

2  2  g 

fc  S2  P  2 

X)  o  a  c 

ts  Cu  X  X 


v-> 

« 

o 

B 

o 

<i> 

> 

6 

o 

B 

p 

1 

c 

'ob 

.5 

u 

ts 

« 

o 

M 

o 

o 

§ 

O 

B 

0) 

Vh 

X 

s 

*0. 

B 

3 

o 

■•4 

</> 

p 

« 

S 

>  OD  00 

§  "Zl  3 


oo  oo  rr 

0o0^c00*^0'x0^p^^c>o 

OoOoOoOo^OOoOO’^O 

II  II  II  II  II  II  II  II  II  II  II  II  II  II  II  II 
<pq<dpQ<dCQ<CQ<;pQ<PQ<Da<PQ 


fS  V  ^ 

1  V)  o  o 


394 


GENERAL  METHODS 


[18] 


[18] 


DYNAMIC  ASSESSMENT  OF  TUMOR  OXYGENATION 


395 


cd 

3 

C/3 


O 

C 

<d 

s 

t/3 

Q 


?S  S 
P!  2 

o  o 


c 

o 

u 

O 

o 


o 

*0 

G 

0)  Q 


c  O 


S.  « 

6-g« 

3  =^1  S 
M  o 

J  CO 


c  2f 
o  c 

U  3 


^  O 


e 

o 

c 

o  .> 

a " 

3  s 

u 

.§ 

ff 

e 

o 

c 

«Q 

JG 

■*-t  Xi 

<c  «a 

oo 

^  -5? 

rt 

43 

pH 

^  P^ 

E 

e:!  s: 

£ 

o  m 

0-1 

rj;  in 

o 

r- 

<s  ^ 

00 

<s 

m  c*-) 

00  S 
P  o 
o  d 


<  PQ  <  P3 


pp 

CQ  PQ 

03  PQ  PQ 

o 

o  o 

o  o  o 

Pn 

Qh  Ph 

tM  Pm  Oh 

Ph 

PM  PM 

PM  PM 

O  v,o 


<  CQ 


^  G 
a>  o 

CA  (A 

O 


4> 


■B 

(U 

o 
G 
(U 
T3 

S  -C 

o' ‘S 

*0  CA 

s  ^ 

3  u 


a  o. 

S  a 

c  .S 

>2  o 

12  o< 
S  <D 
2  <0 
O  43 

a>  o 
o  a 
^  S 

^  C 

.s  O 


•=  o 

r  -  ^ 

»  g  4 
■2.-S  S 

2  p  £ 

g  S  § 

jg  § 

°fe| 

i "  I 

C  o  £ 


fl  r 

Oo 

‘d 

C<3  ro 

•e  c 
S  I 

5? 

S  « 

S  T3 

S  g 

<D  .G 
CL,  O 


Os 

Os 


Os 


✓ 

00 

.  On 


I 


.  ^ 
h  N 


<u 


<D 


I  2  ■  ■& 

■S^-S 

-  P3 

m  ^ 


I 

O 

(5 

S'  ^ 

S'  g 

I 

^  o 

or-  ' 
1^5, 
<§  ro"  2 
.  V)  cc 

n  II 

"o  .S'  ^ 

^5  Cq  'S 


K  C-)  ■ 

1  'o  So 
00 
S  2; 


03 

5  :§ 


G 
O 
u 
O 

,  ^ 
G 

o  ^ 


cL" 
_.  o 
ca  ^ 
G  2 


S  § 
o  cr 
00  (1> 


<  g 

II  .22 


acr;  u 

lA  «  >C] 


^  4: 

CL.  O 

O  tA 

w;  D  S  <( 
1-;  Cl,  T3 

1^^  0  es  « 

§  "c  .iS"  "s 

a  s  3  ^ 

E 

gj  ^  pH  ^ 

J— I  •->  W 


^  S.  R.  Thomas,  L.  C.  Clark,  Jr.  et  al.,  J.  Comput.  Assist.  Tomogr.  10, 1  (1986). 
*  C.  H.  Sotak,  P.  S.  Hees  et  al.,  Magn.  Reson.  Med.  29,  188  (1993). 


396 


GENERAL  METHODS 


[18] 


Ri  is  sensitive  to  temperature,  although  the  response  varies  greatly  be¬ 
tween  PFCs  and  between  individual  resonances  of  each  individual  PFC. 
Over  small  temperature  ranges,  a  linear  correction  to  calibration  curves 
is  appropriate,  but  over  larger  temperature  ranges,  the  response  can  be 
complex,  as  investigated  extensively  by  Shukla  et  a/.’'*”  for  several  PFCs. 
Differential  sensitivity  of  pairs  of  resonances  to  p02  and  temperature 
allowed  Mason  et  to  simultaneously  determine  both  parameters  by 
solving  simultaneous  equations.  However,  generally  it  is  preferable  for  a 
PO2  sensor  to  exhibit  minimal  response  to  temperature,  since  this  is 
not  always  known  precisely  in  vivo  and  temperature  gradients  may  occur 
across  tumors.  As  shown  in  Table  II,  even  a  relatively  small  error  in  tem¬ 
perature  estimate  can  introduce  a  sizable  discrepancy  into  the  apparent 
PO2;  for  example,  the  relative  error  introduced  into  a  p02  determination 
by  a  1°  error  in  temperature  estimate  ranges  from  8  torr/°  for  PFTB''*' 
to  3  torr/°  for  PFOB  (perflubron)’"*^  or  ISCS”’  and  0.1  torr/°  for  HFB,^"^^ 
when  PO2  is  actually  5  torr.  It  must  be  noted  that  error  depends  on  actual 
PO2  and  the  error  varies  with  magnetic  field  and  temperature.  Ri  response 
does  depend  on  magnetic  field,  necessitating  calibration  curves  for  each 
type  of  magnet  system  (e.g.,  1.5,  4.7,  or  7  T).  Thus  comparison  of  PFC 
utility  for  p02  measurements  is  complicated  by  the  field  used  for  specific 
published  investigations,  and  in  Table  II,  we  consider  sensitivity  as 
presented. 

Choice  of  PFC  may  be  governed  by  practical  considerations,  such  as 
cost  and  availability,  since  several  products,  particularly  proprietary  emul¬ 
sions,  may  be  difficult  to  obtain.  HFB  and  15C5  offer  the  immediate  advan¬ 
tage  of  a  high  symmetry  and  a  single  ^®F  NMR  resonance.  This  offers 
maximum  SNR  and  simplifies  imaging,  which  may  otherwise  require  fre¬ 
quency  selective  excitation,  deconvolution,  or  other  NMR  tricks  to  avoid 
chemical  shift  artifacts. 

Route  of  Administration.  The  most  popular  route  for  the  delivery  of 
PFCs  is  as  emulsions  injected  intravenously.  Given  the  extremely  hydro- 
phobic  nature  of  PFCs,  they  do  not  dissolve  in  blood  directly,  but  may  be 
formulated  as  biocompatible  emulsions.  Much  effort  has  been  applied 
to  formulate  stable  homogenous  emulsions,  as  reviewed  elsewhere. 
Following  IV  infusion,  the  emulsion  circulates  in  the  vasculature  with  a 
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typical  half-life  of  12  h,  (depending  on  the  nature  of  the  emulsion)  provid¬ 
ing  substantial  clearance  within  2  days.^"*^’’"^*’  Primary  clearance  is  by 
macrophage  activity,  leading  to  extensive  accumulation  in  the  liver,  spleen, 
and  bone  marrow.^  Indeed,  this  is  a  major  shortcoming  of  IV  delivery, 
since  animals  may  exhibit  extensive  hepatomegaly  or  splenomegaly.^'*®  The 
emulsions  are  not  toxic,  and  other  than  causing  swelling,  appear  not  to 
cause  health  problems.  PFC  clearance  occurs  from  the  liver  with  a  typical 
half-life  of  60  days  for  perfluorotripropylamine  and  3  days  for  perflubron, 
with  primary  clearance  by  migration  to  the  lungs  and  exhalation.*'*® 

Some  investigators  have  examined  p02  of  tissues,  while  PFC  remained 
in  the  blood,  providing  a  vascular  p02.**'*’**®’*^°  Flow  can  generate  artifacts, 
and  correction  algorithms  have  been  proposed.*^*  Many  investigators  have 
measured  p02  in  liver,  spleen,  and  tumors  following  clearance  from  the 
blood,  thus  providing  measurement  of  tissue 

Both  spectroscopic  and  imaging  approaches  have  been  applied  to  tissue 
p02  measurements  depending  on  the  available  SNR.  It  appears  that  uptake 
and  distribution  efficiency  vary  with  tumor  type,  but  in  general,  maximum 
signal  is  detected  from  the  tumor  periphery  corresponding  with  regions  of 
greater  perfusion.*®”  *  Several  reports  have  examined  changes 

in  tumor  p02  in  response  to  acute  interventions  such  as  vasoactive  drugs 
and  hyperoxic  gases.^^’®^’**’®'*^^-*^’’*'’^’**’®  Spectroscopic  time  resolution  has 
ranged  from  seconds  to  minutes,'*’*’**’^  whereas  imaging  often  takes  longer.*^® 
Long  tissue  retention  facilitates  chronic  studies  during  tumor  develop¬ 
ment,  and  progressive  tumor  hypoxiation  has  been  observed  over  extended 
time  periods  of  many  days. *■’'*’*  ^^  Correlated  *®F  and  proton  MRI  suggest 
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that  PFC  does  not  redistribute,  but  remains  associated  with  specific  tissues, 
analogous  to  tree  rings.^^"^’^^^  Thus  in  principle,  a  whole  tumor  can  be  in¬ 
vestigated  by  administering  successive  doses  of  PFC  emulsion  during 
growth. 

PFC  emulsions  may  also  be  administered  intraperitoneally  (IP), 
resulting  in  similar  distribution  to  IV  administration  (unpublished  observa¬ 
tions).  Given  the  volatile  nature  of  many  PFCs,  they  could  be  inhaled,  but 
although  this  is  a  popular  route  for  delivery  of  anesthetics  and  blood  flow 
tracers,  it  does  not  appear  to  have  been  widely  exploited  for  oximetry. 
Nonetheless,  aerosols  have  been  delivered  to  the  lungs  by  inhalation  to 
facilitate  p02  measurements.^^^ 

Two  approaches  have  been  applied  to  circumvent  reticuloendothelial 
uptake.  PFC  has  been  incorporated  in  polyalginate  beads  for  direct  im¬ 
plantation  at  a  site  of  interest.^®"*’*®^  We  favor  direct  IT  injection  of  neat 
PFC,  allowing  any  region  of  interest  in  a  tumor  to  be  interrogated  immedi¬ 
ately.  Use  of  a  fine  needle  ensures  minimal  tissue  damage,  as  described  in 
detail  in  a  later  section.  Others  have  used  direct  injection  of  emulsions  into 
tumors,  but  this  increases  the  volume  considerably,  making  it  more  inva- 
sive.^®  Investigators  have  suggested  that  emulsification  improves  retention 
at  the  site  of  injection.  Direct  injection  of  neat  PFC  also  has  been  used  to 
investigate  retinal  oxygenation'®*’’^®  and  cerebral  oxygenation  in  the 
interstitial  and  ventricular  spaces.’ 

As  described  in  the  following  section  we  favor  direct  intratumoral  injec¬ 
tion  of  neat  HFB  followed  by  echo  planar  imaging  to  generate  p02  maps  in 
tumors. 


FREDOM  (Fluorocarbon  Relaxometry  using  Echo  Planar  Imaging  for 
Dynamic  Ojg^gen  Mapping) 

Recognizing  that  tumors  are  heterogeneous  and  that  p02  may  fluctuate, 
we  developed  a  procedure,  which  allows  repeated  quantitative  maps  of 
regional  p02  to  be  achieved  with  multiple  individual  locations  simulta¬ 
neously  in  6.5  min  with  a  precision  of  1-3  torr,  when  p02  is  in  the 
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Fig.  1.  Hexafluorobenzene  (HFB)  is  a  perfluorocarbon  (PFC)  exhibiting  extensive 
symmetry. 


range  0-10  torr.^®”  We  have  applied  FREDOM  to  diverse  tumor  types  and 
interventions,  as  reviewed  in  a  later  section. 

MRI  is  attractive  because  it  is  readily  available  at  many  institutions.  For 
small  animal  work,  ^®F  NMR  is  widely  available  at  4.7,  7,  and  9.4  T  by 
minor  adaptation  of  routine  instrumentation,  [e.g.,  retuning  proton  radio¬ 
frequency  (RF)  coils].  Within  the  recent  past  ^®F  MRI  is  also  becoming 
available  on  clinical  systems,  facilitating  translation  of  these  techniques  to 
patients.  '^F  NMR  is  particularly  facile  because  there  is  essentially  no 
background  signal  in  tissues  to  interfere  with  measurements,  yet  the  reso¬ 
nance  frequency  and  sensitivity  approach  that  of  proton  NMR.  The  pio¬ 
neering  work  of  Thomas' showed  that  tissue  p02  could  be  imaged  in 
various  organs  based  on  the  '^F  NMR  spin-lattice  relaxation  rate  (Ri)  of 
PFC  reporter  molecules  following  IV  infusion.  Prompted  by  these  studies, 
we  surveyed  a  number  of  PFCs  and  identified  that  HFB  (Fig.  1)  has  many 
virtues  as  a  pOz  reporter.''*'  Symmetry  provides  a  single  narrow  '^F  NMR 
signal,  and  the  spin-lattice  relaxation  rate  is  highly  sensitive  to  changes 
in  pOi,  yet  minimally  responsive  to  temperature  (Fig.  2).''*'  HFB  also  has 
a  long  spin-spin  relaxation  time  (T2),  which  is  particularly  important 
for  imaging  investigations.  From  a  practical  perspective,  HFB  is 
cheap  (<$2/g)  and  readily  available  commercially  in  high  purity  (>99%). 
We  obtain  supplies  from  Lancaster  Synthesis  (Windham,  NH),  though 
many  other  fine  chemical  supply  houses  also  offer  HFB.  We  do  favor 
bottles  over  sealed  ampoules,  since  they  are  easier  to  handle.  HFB  is  well 
eharacterized  in  terms  of  lack  of  toxicity,'®^’'™  exhibiting  no  mutage¬ 
nicity,"'  teratogenicity,  or  fetotoxicity,'’'  and  the  manufacturer’s  material 

Y.  S.  Gorsman  and  T.  A.  Kapitonenko,  Izv.  Estestvennonauchu.  Inst.  Pevinsk.  15,  155 

(1973). 

I.  M.  C.  M.  Rietjens,  A.  Steensma  et  al,  Eur.  J.  Pharmacol.  293,  292  (1995). 

K.  E.  Mortelmans  and  V.  F.  Simmon,  Gov.  Rep.  Announce.  Index  (US)  81,  2555  (1981). 

K.  D.  Courtney  and  J.  E.  Andrews,  J.  Environ.  Sci.  Health  B  19,  83  (1984). 
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A  B 


pOi  (torr)  Temperature  (°) 

Fig.  2.  (A)  The  NMR  spin-lattice  relaxation  rate  (Ri)  of  FIFE  shows  a  linear  response 
to  pOj;  at  37°  and  4.7  T,  Ri  =  0.0835  +  0.001876  x  p02.  (B)  The  ’®F  NMR  Ri  of  HFB  shows 
minimal  response  to  temperature  in  the  physiologic  range;  separate  curves  are  shown  at  four 
different  pO^s:  0,  1,  9.8,  and  21%  O2. 


data  safety  sheet  indicates  LD50  >  25  g/kg  (oral — rat)  and  LC50  95  g/m^/2  h 
(inhalation — mouse). 

HFB  had  been  proposed  as  a  veterinary  anesthetic  and  has  been  used  in 
many  species,  including  ponies,  sheep,  cats,  dogs,  rats,  and  mice,  but  was 
abandoned  because  of  its  high  volatility  (boiling  point  81°)  and  low  flash 
point  (10°). This  presented  unacceptable  dangers  in  veterinary  suites 
for  inhalation  anesthetics.  It  is  not  a  problem  in  our  studies,  where  small 
quantities  of  liquid  (typically,  50  frl)  are  injected  directly  into  the  tumor. 
HFB  requires  no  special  storage,  other  than  a  sealed  bottle  to  prevent 
evaporation.  Melting  point  is  4-6°,  and  density  is  1.62  g/ml. 

Methodology 

Tumor  Preparation.  We  have  applied  this  technique  to  tumors  in  rats 
and  mice,  but  our  methodology  will  focus  on  our  standard  application  to 
rats.  As  for  electrode  techniques,  the  tumor  must  be  accessible  (e.g.,  sub¬ 
cutaneous  in  flank  of  thigh),  though  we  favor  the  pedicle  model,  which  pro¬ 
vides  a  tumor  remote  from  the  body,  analogous  to  an  additional  limb.^^'* 
This  model  is  optimal  for  selective  therapy  such  as  local  hyperthermia,  ir¬ 
radiation,  or  excision.  Noting  the  increasing  interest  in  orthotopic  tumors, 
the  approach  is  also  facile  in  breast  tumors  in  the  mammary  fat  pad  and 
applicable  in  the  prostate  with  a  little  practice. 

L.  W.  Hall,  S.  R.  K.  Jackson  et  al,  in  “Recent  Progress  in  Anaesthesiology  and 

Resuscitation”  (A.  Arias,  R.  Llaurado,  M.  A.  Nalda,  and  J.  N.  Lunn,  eds.),  p.  201.  Excerpta 

Medica,  Oxford,  UK,  1975. 

E.  W.  Hahn,  P.  Peschke  et  al,  Magn.  Reson.  Imaging  11,  1007  (1993). 
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Rats  are  preanesthetized  with  ketamine  hydrochloride  (100  mg/ml)  as 
a  relaxant  and  maintained  under  general  gaseous  anesthesia  with  air 
(1  dm^/min)  and  1.2%  isoflurane  [we  do  note  that  the  appropriate  concen¬ 
tration  of  anesthesia  may  depend  on  strain;  for  example,  Copenhagen  rats 
require  (and  tolerate)  higher  concentrations  than  Fisher  rats].  In  our  latest 
refinement  of  the  procedure,  isoflurane  may  be  stopped  for  a  period  of 
minutes  during  HFB  administration,  since  rats  occasionally  exhibit  respira¬ 
tory  distress,  which  may  be  caused  by  the  anesthetic  properties  of  HFB 
interacting  with  isoflurane.  Within  a  few  minutes,  rats  are  stable  and  rou¬ 
tine  isoflurane  anesthesia  is  maintained.  In  our  earlier  work  we  used 
0.5%  methoxyflurane  in  33%  oxygen  with  66%  but  isoflurane 

appears  to  be  a  less  stressful  anesthetic. 

HFB  is  deoxygenated  by  bubbling  nitrogen  for  5  min  before  use.  We 
have  previously  shown  in  hypoxic  tumor  biopsies  that  use  of  aerated 
HFB  may  introduce  a  systematic  apparent  elevation  in  p02  (~2-3  torr). 
HFB  is  injected  directly  into  tumors  using  a  gas  tight  syringe  with  a 
custom-made  fine  sharp  needle  (32G  #7803-04;  Hamilton,  Reno,  NV). 
Generally,  HFB  is  administered  along  three  to  five  tracks  in  the  form  of 
a  fan  in  a  single  central  plane  of  the  tumor  coronal  to  the  rat’s  body 
(Fig.  3A).  The  needle  is  inserted  manually  to  penetrate  across  the  whole 
tumor  and  withdrawn  ~1  mm  to  reduce  pressure,  and  3  pi  HFB  is 
deposited.  The  needle  is  repeatedly  withdrawn  a  further  2-3  mm,  and 
additional  HFB  is  deposited.  Typically,  HFB  is  deliberately  deposited  at 
about  16  individual  locations  per  tumor,  in  both  the  central  and  pe¬ 
ripheral  regions  of  the  tumors,  to  ensure  that  the  interrogated  regions  are 
representative  of  the  whole  tumor. 

The  animal  is  placed  on  its  side  in  a  cradle  with  a  thermal  blanket  to 
maintain  body  temperature.  A  fiber-optic  probe  is  inserted  rectally  to 
monitor  core  temperature.  Temperature  measurement  is  optional,  since 
the  Rj  is  essentially  invariant  with  temperature  and  does  not  need  correc¬ 
tion  for  p02  estimates.  Of  course,  temperature  regulation  and  measure¬ 
ment  is  important  to  ensure  stable  tumor  physiology.  An  MR-compatible 
pulse  oximeter  equipped  with  “rat”  software  (8600V;  Nonin,  Inc.,  Ply¬ 
mouth,  MN)  may  be  applied  to  a  hind  foot  to  monitor  arterial  oxygenation 
(Sa02)  and  heart  rate  (optional,  but  provides  additional  useful  data 
regarding  animal  health  and  physiology). 

Most  of  our  MR  experiments  were  performed  using  an  Omega  CSI  4.7 
horizontal  bore  magnet  system  with  actively  shielded  gradients  (GE 
systems,  acquired  by  Bruker  Instrument,  Inc.,  Fremont,  CA).  Recently, 
our  MR  system  has  been  upgraded  to  a  Varian  Unity  INOVA,  providing 

R.  P.  Mason,  A.  Constantinescu  et  al.,  Radiat.  Res.  152,  239  (1999). 
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Fig.  3.  (A)  Recommended  pattern  of  IT  injection.  Several  tracks  in  the  form  of  a  fan  in  a 
single  plane.  (B)  Four  contiguous  MR  images  (5  mm  thick)  of  a  representative  Dunning 
prostate  R3327-AT1  tumor  (volume  =  1.6  cm^).  Top:  MRI;  middle:  corresponding  '®F; 

bottom:  overlay  to  show  interrogated  regions.  B  marks  point  of  attachment  of  the  tumor 
to  back  of  rat. 
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- —  MBEST 

RF  pulses 
G,(Slice) 


Fig.  4.  MRI  sequence  used  to  achieve  relaxation  curves. 


enhanced  data  acquisition  and  processing  software  and  stronger  imaging 
gradients.  A  tunable  (^H/^^F)  MR  coil,  2  or  3  cm  in  diameter  matched  to 
the  tumor  size  (constructed  from  a  cylindrical  copper  tube  about  2  cm  deep 
and  acting  as  a  single-turn  solenoid),  is  placed  around  the  tumor-bearing 
pedicle.  Shimming  is  performed  on  the  proton  water  signal  to  a  typical  line- 
width  of  60-100  Hz.  Proton  images  are  obtained  for  anatomic  reference 
using  a  three-dimensional  (3D)  spin-echo  sequence.  The  coil  is  then  re¬ 
tuned  in  place  to  188.27  MHz,  and  corresponding  MR  images  are 
obtained.  Overlaying  the  ^^F  MR  images  on  the  corresponding  proton 
images  reveals  the  distribution  of  HFB  (Fig.  3B).  Typically,  ^®F  NMR 
signal  is  obtained  from  6-10%  of  the  total  tumor  voxels.^^®'^^^'^^^ 

p02  maps  are  obtained  using  our  standard  MR  oximetry  proto- 
coi  116,160.176  jjjjg  applies  pulse  burst  saturation  recovery  (PBSR)  echo 
planar  imaging  (EPI)  relaxometry  using  the  ARDVARC  (alternated  relax¬ 
ation  delays  with  variable  acquisitions  to  reduce  clearance  effects)  acquisi¬ 
tion  protocol  to  map  the  tumors  (Fig.  EPI  uses  a  single  spin-echo 

with  “blipped”  phase  encoding  (modulus-blipped  echo-planar  single-pulse 
technique  [MBEST]),  although  other  EPI  sequences  should  be  equally  ap¬ 
plicable.  We  chose  the  PBSR  approach  to  Ti  relaxation  measurements  for 
historical  reasons;  our  earliest  work  used  ^^F  NMR  spectroscopy  of  tumors 
with  excitation  and  detection  based  on  surface  coils.^^^  The  PBSR  ap¬ 
proach  is  suitable  for  use  with  the  nonuniform  excitation  typical  of  surface 


B.  R.  Barker,  R.  P.  Mason  et  al,  Magn.  Resort.  Imaging  11,  1165  (1993). 
D.  Le,  R.  P.  Mason  et  al,  Magn.  Reson.  Imaging  15,  971  (1997). 
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Fig.  5.  Typical  relaxation  curve  showing  order  of  data  acquisition,  y  =  A  [1  -  (1  +  W) 
exp(-Ri  X  t)]  and  Tj  =  10.34  ±  0.16  s,  A  =  99.5,  W  =  0.04,  and  hence  p02  =  7.0  ±  0.7  torr. 


coils  as  presented  by  Evelhoch  and  Ackerman.^’®  For  EPI,  we  do  apply 
carefully  calibrated  nil  pulses  to  ensure  accurate  refocusing  of  images, 
but  we  continue  to  favor  the  PBSR  approach,  since  it  provides  particularly 
rapid  estimation  of  Ti  by  avoiding  the  need  for  extended  relaxation  recov¬ 
ery  times  between  acquisitions.  Saturation  is  achieved  by  a  series  of  20  non- 
spatially  selective  ir/2  pulses  with  50-ms  spacing,  which  is  sometimes  called 
a  “Comb”  format. 

Typical  FREDOM  parameters  use  32  x  32  data  points  across  a  field  of 
view  of  40  X  40  mm,  providing  1.25  mm  in-plane  resolution.  Recently,  we 
have  applied  a  slice  selection  gradient,  providing  5-mm-thick  slices,  and 
hence,  8-/xl  voxels.  However,  the  deposition  of  HFB  is  designed  to  occur 
in  a  plane,  and  thus  the  discrete  distribution  of  the  reporter  molecule  itself 
can  be  used  to  define  the  slice.  We  apply  a  PBSR  preparation  sequence  be¬ 
cause  it  is  ideally  suited  for  measuring  diverse,  long  Ti  values.  Unlike  more 
traditional  inversion  recovery  sequences,  there  is  no  need  to  wait  >5  Ti  be¬ 
tween  successive  images.  As  shown  in  Fig.  5,  we  use  14  delays  in  the  order 
90  s,  200  ms,  60  s,  400  ms,  40  s,  600  ms,  20  s,  800  ms,  16  s,  1  s,  8  s,  1.5  s,  4  s, 
and  2  s,  selected  to  cover  the  whole  range  of  TjS  (viz  p02  values).  Many 
papers  have  been  published  on  optimizing  relaxation  curves  and  choosing 
parameters  to  enhance  precision.^’®^’*^^  Our  experience  suggests  that  our 


J.  L.  Evelhoch  and  J.  J.  H.  Ackerman,  /.  Magn.  Resort.  53,  52  (1983). 
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parameters  are  appropriate  for  interrogating  a  broad  range  of  pOa  values, 
though  they  have  not  been  rigorously  or  theoretically  optimized.  More 
delays  could  improve  curve  fitting,  but  would  increase  experiment  time. 
We  believe  fewer  delays  would  degrade  quality  of  the  curve  fitting. 

Traditional  Tj  measurement  sequences  acquire  data  with  delays  in 
monotonic  order,  whereas  we  alternate  longer  and  shorter  delays  to  mini¬ 
mize  any  systematic  errors,  which  would  be  introduced  if  the  signal  ampli¬ 
tude  varies  during  the  measurement.  We  have  found  that  HFB  clears  from 
tissue  with  a  typical  half-life  of  600  min,^®^  which  would  introduce  errors 
into  the  amplitude.  Although  the  total  acquisition  time  for  a  Ti  map  is 
6.5  min,  we  reduce  the  time  between  first  and  last  acquisitions  further 
by  applying  the  longest  delay  first.  Our  experience  shows  that  it  is  import¬ 
ant  to  measure  at  least  two  data  points  on  the  plateau  of  the  curve  (i.e., 
>5  X  Ti,  and  we  choose  60  and  90  s,  respectively.  The  data  points  with 
longer  recovery  times  have  greater  SNR,  and  we  find  that  a  minimum 
SNR  =  10  is  required  to  produce  satisfactory  Ti  curves  and  p02  estimates. 
The  variation  in  amplitude  between  longest  and  shortest  delays  approaches 
100-fold.  The  poorest  SNR  data  points  (short  delay;  r)  could  compromise 
the  quality  of  the  Ti  curves,  and  thus  we  obtain  multiple  acquisitions  at 
these  times  to  provide  enhanced  SNR  by  signal  averaging.  Because  the  r 
values  are  short,  it  adds  little  to  the  overall  experimental  time.  We  use 
NA  =  12  for  r  <  1  s;  NA  =  8  for  r  =  1  s;  NA  =  4  for  r  =  1.5,  2,  and  4  s; 
NA  =  2  for  T  =  8  s;  and  16  s  and  NA  =  1  for  the  longer  delays.  Signal  am¬ 
plitudes  are  corrected  for  the  additional  acquisitions,  and  a  2-d  Fourier 
transform  (FT)  is  applied  to  each  image.  Curve  fitting  is  then  applied  on 
a  voxel-by-voxel  basis  to  each  image  set.  Data  quality  could  probably  be 
enhanced  by  application  of  apodization  and  filtering  functions,  such  as 
Fermi  and  Hanning  filters,  though  we  have  not  yet  implemented  these 
approaches.  The  acquisition  protocol  is  not  part  of  the  standard  software 
supplied  by  the  manufacturer,  but  we  can  assist  interested  investigators. 

Data  are  transferred  to  a  personal  computer  (PC)  for  further  analysis 
using  a  program  created  in  our  laboratory  (Timap  [Fig.  6],  PASCAL  by 
Dennis  Le).  Timap  recognizes  issues  in  data  analysis  and  provides  filters, 
clustering  algorithms,  and  temporospatial  correlation  to  assist  in  effective 
data  reduction.  At  top  right  of  Fig.  6  is  a  map,  which  may  be  toggled  to  dis¬ 
play  signal  amplitude,  Ti,  Ti  error  (Ti  err),  or  p02  (Fig.  7).  Blank  pixels  in¬ 
dicate  that  no  curve  fit  could  be  achieved.  Gray  pixels  indicate  a  successful 
curve  fit  and  potential  p02  value,  but  with  errors  beyond  the  specified 
range.  In  this  example,  thresholds  were  set  as  Ti  err  <  2.5  s  and  the  ratio 
Ti  err/Ti  <  50%.  Tightening  these  thresholds  can  produce  higher-quality 
data,  but  eliminates  more  data.  Colored  pixels  report  measurements 
satisfying  the  threshold  criteria,  and  these  are  tabulated  in  the  form  of  a 
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Fig.  6.  Example  Timap  used  for  data  filtering,  clustering,  and  analysis. 
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histogram  (center),  which  can  be  plotted  with  selected  data  bins  and  ranges 
of  interest.  At  bottom  right  of  Fig.  6  are  data  for  each  pixel,  including  co¬ 
ordinate  and  fitted  parameters.  In  the  box  on  the  left  are  shown  the  number 
of  fitted  data  points,  number  of  accepted  data  points,  and  statistics  such  as 
mean  Ti,  pOa,  and  corresponding  amplitude-weighted  values. 

pOa  is  calculated  on  the  basis  of  a  curve  fit  to  exponential  data.  Any 
curve  fitting  is  associated  with  uncertainty,  and  indeed  our  procedure  deter¬ 
mines  Ti  err  values.  Although  the  concept  of  negative  pOa  values  may  appear 
impossible,  it  is  perfectly  reasonable,  provided  that  error  bars  and  uncer¬ 
tainties  are  considered.  Thus  a  value  of  -2  torr  ±  3  torr  legitimately  indi¬ 
cates  an  actual  pOa  very  close  to  zero.  Some  oximetry  approaches  ignore 
negative  values  or  bin  them  all  as  zero.  We  accept  any  value  providing  that 
the  error  estimated  for  Ti  is  within  a  specified  range. 

The  most  powerful  aspect  of  FREDOM  is  the  ability  to  follow  the  fate 
of  individual  voxels,  and  thus  we  usually  acquire  at  least  three  baseline  pOa 
maps,  followed  by  further  maps  accompanying  interventions,  such  as  hy- 
peroxic  gas  breathing.  Even  under  baseline  conditions,  fluctuations  in  Tj 
are  apparent.  These  may  arise  from  uncertainty  in  Ti,  which  may  be 
reflected  in  Ti  err  or  transient  fluctuations  in  pOa. 

Ti  map  allows  us  to  investigate  sequential  maps  and  select  only  those 
pixels  that  consistently  show  small  Ti  en-s,  for  example,  white  regions  of 
interest  (ROI).  In  subsequent  maps  we  may  include  only  these  regions, 
which  satisfy  the  inclusion  criteria  for  every  map.  This  allows  us  to  follow 
population  dynamics  based  on  specific  tumor  regions,  avoiding  potential 
anomalies  due  to  varying  numbers  of  “good  quality”  pixels  in  sequential 
images.  Typically,  50^150  pixels  provide  high-quality  data  in  any  given 
map,  but  generally  20-80  pixels  may  be  followed  for  a  series  of  23  maps 
associated  with  interventions. 

The  ROI  tool  also  facilitates  clustering,  for  example,  selection  of  only 
those  tumor  regions,  which  consistently  have  pOa  <  10  torr  throughout 
the  baseline  period  (these  may  be  considered  as  chronically  hypoxic,  as 
opposed  to  regions,  which  fluctuate  and  may  only  be  <10  torr  in  some 


Fig.  7.  Maps  showing  (A)  '®F  NMR  signal  intensity,  (B)  spin-lattice  relaxation  times  (Tj), 
(C)  estimated  errors  in  Tj,  and  (D)  p02  values.  Black  squares  show  regions  where  signal  to 
noise  or  data  quality  was  so  poor  as  to  not  provide  a  T i  curve  fit.  Gray  voxels  show  regions  that 
provided  a  curve  fit,  but  where  the  uncertainty  in  the  data  exceeded  the  threshold  criteria  (i.e., 
'T'l  err  ^  2*5  s  or  Ti  errTi  >  50%).  All  colored  voxels  (n  =  100)  provided  a  curve  fit  within  the 
acceptance  criteria  in  this  map,  but  only  those  56  voxels  within  the  white  regions  of  interest 
(ROI)  provided  consistently  high-quality  data  throughout  the  sequence  of  eight  maps  acquired 
during  baseline  and  with  oxygen  breathing.  Data  were  obtained  using  the  FREDOM  approach 
from  a  representative  Dunning  prostate  R3327-AT1  tumor  (2  cm^). 
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baseline  maps).  The  ability  to  follow  groups  of  pixels  with  particular  base¬ 
line  characteristics  has  revealed  heterogeneity  in  response  to  many  inter¬ 
ventions — often  those  regions  initially  well  oxygenated  show  rapid  and 
large  response  to  hyperoxic  gas  breathing,  whereas  those  that  are  initially 
poorly  oxygenated  show  little  and  sluggish  response  (e.g.,  Fig.  8).  This  ap¬ 
proach  can  also  reveal  phenomena  such  as  the  “steal”  effect,  whereby  ini¬ 
tially  well-oxygenated  regions  decrease  in  p02,  while  others  increase, 
which  might  appear  as  “no  change,”  when  histogram-based  population 
statistics  are  used. 

Figure  9  shows  that  the  data  quality  is  strongly  related  to  signal  ampli¬ 
tude.  Although  these  data  represent  56  voxels  from  a  single  tumor,  we  have 
previously  shown  similar  data  for  the  alternate  PFC  perflubron.^’  ®  Given  the 
high  solubility  of  O2  in  PFCs,  there  could  be  a  concern  that  PFCs  act  as 
reservoirs,  perturbing  local  p02.  Fig.  9B  shows  that  there  is  no  correlation 
(r^  <  0.05)  between  signal  amplitude  (viz.  HFB  concentration)  and  p02. 

As  with  any  measurement,  sampling  is  a  critical  issue.  FREDOM  is 
analogous  to  the  Eppendorf  Histograph  in  that  it  samples  multiple  loca¬ 
tions,  which  appear  to  reflect  interstitial  p02.  Histograph  data  suggest  that 
a  minimum  of  100-140  data  points^®^  along  five  tracks^®^  are  required  to  ac¬ 
curately  represent  the  p02  distribution  of  a  tumor,  though  such  criteria  do 
depend  on  tumor  size  and  heterogeneity.  We  applied  Monte  Carlo  simula¬ 
tion  to  assess  the  data  requirements  for  FREDOM.  A  data  set  was 
selected,  which  provided  120  high-quality  data  points,  and  these  data  points 
were  accessed  in  random  order  with  continuous  calculation  of  mean  p02  to 
assess  the  asymptotic  trend  lines  toward  the  actual  mean  p02  (Fig.  10).  It 
appears  that  about  50  data  points  are  required  to  well  represent  the  p02 
distribution  of  the  tumor,  which  is  generally  achievable  in  any  individual 
PO2  map.  However,  the  unique  ability  to  observe  dynamic  changes  in 
PO2  at  multiple  locations  simultaneously  is  the  greatest  strength  of  FRE¬ 
DOM.  Detection  of  heterogeneous  responses  is  useful  even  if  fewer  data 
points  are  examined,  since  each  location  serves  as  its  own  control. 

Validation  of  Measurements.  The  spin-lattice  relaxation  rate  (Ri)  is  de¬ 
termined  for  each  voxel  using  a  three-parameter  fit  of  signal  intensities  to 

y;  =  A[1  -  (1  -f-  W)  exp(-Ri  x  r)]  (4) 

using  the  Levenberg-Marquardt  least  squares  fitting  protocol.^  ’*’  Typically, 
~100-300  voxels  provide  an  Ri  fit  and  potential  p02  value.  Because  noise 
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Fig.  9.  (A)  Relationship  of  Rj  with  signal  amplitude  showing  strong  correlation  (r^  > 
0.69).  (B)  Relationship  of  p02  and  signal  amplitude  showing  lack  of  correlation  (r^  <  0.05). 


Fig.  8.  PO2  maps  obtained  using  the  FREDOM  approach  from  the  tumor  shown  in  Fig.  3. 
(A)  Under  baseline  conditions,  those  voxels  (n  =  56)  within  the  white  region  provided 
consistently  reliable  data  during  repeated  measurements.  (B)  Baseline  map  (breathing  air: 
FO2  =  21%);  mean  p02  =  7.2  ±  2.6  (SE)  torr,  median  p02  =  1.3  torr  (range  -7-88  torr). 
(C)  Breathing  oxygen  (FO2  =  100%):  fourth  map  obtained  24-32  min  after  switching  from  air: 
mean  p02  =  47.2  ±  10.7  torr  {p  <  0.0001  compared  with  baseline),  median  p02  =  8.3  torr 
(range  -3-204  torr).  (D)  Breathing  carbogen  (FO2  =  95%):  fifth  map  after  switching  to 
carbogen:  mean  p02  =  43.1  ±  9.2  torr  (p  <  0.0001  compared  with  baseline),  median  p02  =  5.6 
torr  (range  -10-216  torr). 
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Fig.  10.  Monte  Carlo  simulation  of  tumor  oxygenation.  Asymptotic  behavior  shows  that 
tumor  may  be  characterized  by  about  50  data  points.  For  this  large  ATI  tumor  in  a  rat 
breathing  oxygen,  all  estimates  converge  on  the  mean  p02  =  31.7  +  1.0  (SE)  torr. 


itself  may  give  an  apparent  relaxation  curve  (Ri)  fit,  data  are  selected 
within  a  region  of  interest,  having  Ti  err  <  2.5  s  and  the  ratio  Tj  err/Ti  < 
50%.  With  respect  to  respiratory  interventions,  only  those  voxels  that  pro¬ 
vided  consistently  reliable  data  throughout  the  measurements  are  included 
for  further  analysis.  At  37°  and  4.7  T 

p02(torr)  =  [Ri(s-i)  -  0.0835]/0.001876  (5) 

Equation  (5)  provides  p02  in  units  of  torr.  The  literature  can  be  compli¬ 
cated  by  use  of  various  units,  and  it  may  be  instructive  to  provide  conver¬ 
sion  factors.  We  favor  torr  (1  torr  =  1  mmHg),  since  radiobiologic  hypoxia 
develops  in  the  range  0-15  torr  and  this  is  the  traditional  unit  favored 
by  radiation  biologists  and  oncologists.  Further,  760  torr  =  1  standard  at¬ 
mosphere  (atm.),  and  gases  are  often  quoted  in  %atm.  For  SI  units,  101, 
325  Pa  (or  N/m^)  =  1.01325  x  10^  dynes/cm^  =  1  atm.  Some  investigators 
quote  oxygen  concentrations  in  fiM,  often  assuming  that  the  solubility  of 
oxygen  in  water  at  37°  is  1.35  iiM HoirP  Quoting  concentrations  can  be 
confusing  because  the  solubility  of  oxygen  is  highly  variable  in  solutions. 
FREDOM  calibration  is  based  on  p02  values,  and  as  shown  by  Eqs.  1-3, 
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these  are  directly  related  to  the  solubility  of  oxygen  in  HFB.  However, 
partitioning  of  oxygen  between  the  aqueous  and  PFC  phases  depends  on 
PO2,  not  concentrations.  Thus  the  p02  determined  by  FREDOM  will  ac¬ 
curately  reflect  the  ambient  tissue  p02  and  if  the  solubility  of  oxygen  can 
be  estimated  in  the  particular  milieu,  then  [O2]  could  be  calculated. 

Both  systematic  and  random  errors  may  interfere.  Random  errors  may 
be  diminished  by  performing  multiple  repeat  measurements  (provided  that 
the  biologic  system  is  stable).  Systematic  errors  are  more  complex  and 
could  arise  inter  alia  from  erroneous  calibrations,  inappropriate  curve  fit¬ 
ting,  and  temperature  changes.  Appropriate  curve  fitting  may  be  the 
greatest  problem  in  relaxation  analysis.  Provided  signal/noise  >  10  for  the 
most  intense  signals,  we  generally  obtain  excellent  curve  fits.  Each  voxel 
(or  ROI)  will  comprise  HFB  at  a  range  of  p02  values,  creating  a  multiex¬ 
ponential  curve.  However,  modeling  shows  that  the  fit  provides  an  “aver¬ 
age”  value.  For  population  data,  we  sometimes  determine  error-weighted 
means  [S(x/it^)]/[S(1/(t2)]  in  order  to  exploit  as  much  available  data  as 
possible. 

In  terms  of  absolute  p02  values,  the  most  critical  aspect  is  effective  cali¬ 
bration  curves.  Over  the  years  we  have  achieved  p02-dependent  ^^F  NMR 
relaxation  curves  for  many  and  encountered  potential 

pitfalls.  The  calibration  curve  we  recommend  at  4.7  T  and  37°  is  given  by 
Eq.  (5)  which  was  originally  presented  by  Hunjan  et  al}^°  Briefly,  125  ^1 
HFB  was  added  to  each  of  four  gas-tight  NMR  tubes  together  with 
0.5  ml  water  and  saturated  at  37°  by  bubbling  with  carbon  dioxide,  1% 
O2  (balance  N2),  9.8%  O2  (balance  N2),  or  air,  respectively.  Tubes  were 
sealed,  and  the  phantom  was  maintained  at  37°  in  a  water  bath  within  a  coil 
in  the  magnet.  FREDOM  was  applied  using  the  parameters  described 
earlier,  and  the  spin-lattice  relaxation  rates  were  estimated  on  a  voxel- 
by-voxel  basis  using  three -parameter  fit.  Equation  (5)  was  established  using 
linear  regression  analysis  of  amplitude  squared  weighted  mean  values  for 
each  gas.^^^ 

Although  the  relationship  Rj  =  f  (PO2)  is  theoretically  expected  to  be 
linear  and  empirically  found  to  be  so,  we  believe  that  it  is  important  to  use 
calibration  gases  in  the  range  of  physiologic  p02.  We  recommend  purchase 
of  rigorously  calibrated  gases.  We  bubble  gases  for  30  min  and  use  gas-tight 
Wilmad  NMR  tubes,  which  may  be  sealed  with  ground  glass  joints. 
Samples  are  used  within  hours  of  saturation.  As  desired,  multiple  tubes 
may  be  prepared  at  given  p02,  but  since  a  linear  relationship  is  expected, 
it  may  be  equally  appropriate  to  use  additional  gases.  Ti  may  be  deter¬ 
mined  multiple  times  for  each  sample,  and  we  recommend  using  the  pulse 
sequence  to  be  applied  for  in  vivo  investigations.  In  1997,  we  published  a 
calibration  curve 
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Ri(s-i)  =  0.074(±0.003)  +  0.012(±0.0002)p02(%atm.)  (6) 

for  HFB  at  37°  and  4.7  Considering  temperature  sensitivity  gave 

Ri(s-i)  =  0.77(±0.03)  -  (0.00009  ±  0.001)TfC) 

+  (0.018  ±  0.003)P(%atm.)  -(0.00017  ±  0.00008)TP(°C%  atm.) 

(V) 

More  recently,  we  achieved  Eq.  (5)  as  a  new  calibration  curve  at  37°  and 
4.7  T  using  the  ARDVARC  approach.’^”  This  calibration  curve  is  based 
on  higher-quality  data  and  provides  superior  p02  estimates.  In  particular, 
when  tissue  is  expected  to  be  hypoxic  (excised  tissue),  we  find  fewer 
apparently  negative  p02  values. 

A  major  strength  of  the  FREDOM  approach  is  that  calibration  curves 
remain  valid  between  samples  and  across  experimental  platforms.  Calibra¬ 
tion  curves  obtained  in  vitro  are  valid  in  vivo.  Thus  we  believe  other  inves¬ 
tigators  can  apply  Eq.  (5)  without  the  need  to  establish  their  own 
calibrations,  provided  that  studies  are  undertaken  at  about  4.7  T  and  in 
the  range  30-45°.  As  discussed  in  an  earlier  section  and  Table  II,  a  1°  error 
in  temperature  estimates  will  only  introduce  a  0.13-torr  error  in  p02  esti¬ 
mate,  when  p02  is  about  5  torr.  We  also  note  that  the  calibration  curves 
for  HFB  are  relatively  insensitive  to  magnetic  field.’ 

We  have  previously  investigated  microbiodistribution  of  HFB  based  on 
Oil  Red  O  stain,  which  indicated  that  HFB  occurs  as  microscopic  droplets 
(1-20  ^m)  widely  distributed  across  tumor  tissue.’’^  Significantly,  there 
was  no  evidence  for  formation  of  films,  which  could  act  as  conduction  con¬ 
duits,  causing  oxygen  equilibration.  Occasionally,  an  animal  (viz.  tumor) 
will  move  slightly  during  a  long  time  series  of  measurements.  In  this  case 
the  requirement  of  consistently  high-quality  curve  fits  throughout  the  data 
set  for  individual  voxels  fails  because  the  tissues  have  moved  relative  to  the 
voxel  grid.  Such  motion  is  immediately  apparent  when  examining  the 
images.  There  is  a  choice  of  eliminating  such  a  data  set  or  relaxing  the  ac¬ 
ceptance  criteria  and  examining  p02  population  distribution  without 
spatial  continuity.  It  is  important  to  recognize  that  all  methods  include  a 
degree  of  sampling.  In  some  approaches,  this  involves  the  selective  place¬ 
ment  and  tracking  of  electrodes;  in  others,  the  choice  of  locations  for 
biopsy  and  the  number  of  microscopic  fields  of  view.  Sampling  can  be 
avoided  by  obtaining  global  measurements,  as  commonly  acquired  using 
near-infrared  approaches,  but  this  may  itself  mask  the  fundamental  tumor 
heterogeneity’®’  (see  Chapter  17  in  this  volume). 

We  have  previously  shown  that  HFB  shows  little  macroscopic  redistri¬ 
bution  over  a  period  of  hours.’^'’  It  does  clear  from  the  tumors  with  a 
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typical  half-life  of  HFB  of  about  600  min,  though  some  tumors  show  essen¬ 
tially  no  detectable  clearance  over  a  period  of  6  Clearance  of  HFB 

precludes  long-term  studies  of  chronic  oxygenation,  unless  further  doses  of 
HFB  are  administered. 

Comparison  of  pOa  distributions  using  FREDOM  or  the  Eppendorf 
Histograph  has  shown  close  similarity  in  both  small  and  large  tumors.^^^ 
Dynamic  studies  in  several  tumor  types  have  shown  equivalent  behav¬ 
ior  when  assessed  using  polarographic  oxygen  electrodes  or  OxyLite  or 
FOXY  optical  probes.^^’'*^’^^''  Relative  hypoxia  has  been  compared  with 
the  histologic  reporter  pimonidazole,  revealing  similar  trends  across  tumor 
types.^^^ 

Applications  of  FREDOM 

FREDOM  has  been  applied  to  investigations  of  diverse  tumor  types 
(syngeneic  rat  prostate  and  breast  tumors  and  xenograft  human  lym¬ 
phomas)  with  respect  to  growth  and  acute  interventions.^^’^^’’^'*"^^^’^®'^  Per¬ 
haps  the  greatest  strength  is  the  ability  to  investigate  regional  dynamic 
changes  in  pOj  accompanying  acute  interventions.  Figure  8  shows  changes 
in  regional  pOa  in  an  ATI  tumor  accompanying  hyperoxic  gas  breathing. 
As  published  previously the  areas  that  were  initially  better  oxy¬ 
genated  responded  with  elevated  pOa,  whereas  those  relatively  hypoxic 
regions  showed  little  response.  Data  also  may  be  presented  as  histograms 
(Fig.  11),  revealing  significant  differences  between  mean  and  median  p02 
and  hypoxic  fractions  between  small  and  large  tumors,  and  between  the 
slow-  and  fast-growing  sublines  H  and  ATI  of  the  Dunning  prostate 
R3327.  Figure  12  shows  variation  in  global  mean  pOa  accompanying  re¬ 
spiratory  challenge  with  oxygen  or  carbogen  and  return  to  baseline 
following  the  intervention.  Because  individual  tumor  regions  may  be  ob¬ 
served,  local  response  may  be  compared  at  identical  locations,  efficiently 
comparing  the  efficacy  of  interventions  (Fig.  13).  Figure  14  shows  the  dif¬ 
ferential  behavior  of  regions  in  undifferentiated  Dunning  prostate  R3327- 
ATl  tumors  versus  highly  differentiated  H  tumors.  In  each  tumor  type 
there  are  both  well  and  poorly  oxygenated  regions  under  baseline  condi¬ 
tions  (see  Fig.  11).  In  response  to  hyperoxic  gas  breathing,  well-oxygenated 
regions  in  both  tumor  types  show  a  rapid  and  significant  elevation  in  p02, 
which  is  reversible  when  inhaled  gas  is  returned  to  air  (Fig.  14).  Poorly  oxy¬ 
genated  regions  in  H  tumors  respond  slowly,  but  ultimately  rise  above  the 
range  of  radiobiologic  hypoxia,  whereas  corresponding  regions  in  ATI 
tumors  do  not.  We  also  have  investigated  vasoactive  drugs'^^  and  vascular 
targeting  agents,^^  and  the  short-term  changes  in  p02  following  irradiation 
also  have  been  examined.^**'*  Perhaps  the  most  significant  results  to  date 
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Fio.  11,  Histograms  of  tumor  oxygenation  pooled  for  six  H  (A)  and  seven  ATI  (B)  small 
Dunning  prostate  R3327  rat  tumors.  The  H  tumors  (219  voxels)  showed  significantly  higher 
baseline  p02  (mean  =  31  ±  2  torr,  median  =  27  torr)  than  the  size-matched  ATI  tumors  (338 
voxels;  mean  =  14  ±  1  torr,  median  =  11  torr;  p  <  0.0001).  With  respect  to  oxygen  challenge 
(bottom),  mean  p02  in  both  groups  increased  significantly  (mean  =  122  ±  7  in  the  H  versus, 
mean  =  58  ±  4  in  the  ATI;  p  <  0.0001).  Arrows  indicate  mean  (x)  and  median  (m)  pOa, 
respectively. 


show  that  p02  measurements  and  detection  of  changes  in  p02  accompany¬ 
ing  interventions  correlate  with  the  efficacy  of  tumor  irradiation.^^® 


Future 

Ultimately,  the  value  of  a  technique  depends  on  its  robustness,  ease  of 
use,  and  widespread  implementation.  To  date,  few  laboratories  had 
adopted  the  FREDOM  approach  because  efficient  investigation  of  HFB 
requires  an  unusual  NMR  pulse  sequence.  With  the  recent  upgrade  of 
our  own  instrumentation  to  the  Varian  Unity  INOVA,  the  software  is 
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Fig.  12.  Dynamic  pOa  (mean  ±  SE)  obtained  from  sequential  maps  of  the  ATI  tumor 
shown  in  Fig.  8  with  respect  to  respiratory  challenge.  *,  p  <  0.001;  **,  p  <  0.0001  versus 
baseline. 


Fig.  13.  Correlation  between  maximum  p02  detected  in  each  of  56  voxels  from  the  ATI 
tumor  shown  in  Fig.  8,  when  the  rat  breathed  oxygen  versus  carbogen  (r^  >  0.85). 


now  available  on  this  popular  platform,  facilitating  ready  implementation 
elsewhere.  In  terms  of  research  applications,  it  is  known  that  tumor  tissue 
p02  varies  rapidly  in  response  to  many  acute  interventions,  ranging  from 
irradiation  to  photodynamic  therapy,  and  various  chemotherapies.  We 
foresee  FREDOM  as  a  valuable  tool  for  assessing  the  dynamic  time  course 
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Fig.  14.  Variation  in  pOa  with  respiratory  challenge  for  individual  regions  chosen  as 
initially  well  oxygenated  (pOa  >  10  torr)  or  hypoxic  (pOj  <  5  torr)  from  an  ATI  (dotted  lines) 
and  H  (solid  lines)  tumor,  respectively.  All  the  well-oxygenated  regions  in  the  two  sublines 
increased  significantly  in  pOa  in  response  to  oxygen  breathing  (p  <  0.01).  The  hypoxic  regions 
from  the  H  tumor  increased,  whereas  those  in  the  ATI  tumor  did  not. 


of  such  interventions  to  provide  clear  insight  into  the  mode  of  action  of 
therapeutic  approaches  and  aid  in  the  high-throughput  screening  of  new 
drugs,  such  as  vascular  targeting  and  antiangiogenic  agents. 
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TUMOR  PHYSIOLOGIC  RESPONSE  TO  COMBRETASTATIN  A4  PHOSPHATE 

ASSESSED  BY  MRI 


Dawen  Zhao,  M.D.,  Ph.D.,  Lan  Jiang,  M.S.,  Eric  W.  Hahn,  Ph.D.,  and  Ralph  P.  Mason,  Ph.D. 

Department  of  Radiology,  University  of  Texas  Southwestern  Medical  Center,  Dallas,  TX 

Purpose;  To  evaluate  the  effect  of  the  vascular  targeting  agent,  combretastatin  A4  phosphate,  on  tumor 
oxygenation  compared  with  vascular  perfusion/permeability. 

Methods  and  Materials;  MRI  oximetry  and  dynamic  contrast-enhanced  (DCE)-MRI  were  used  to  monitor 
tumor  oxygenation  and  perfusion/permeability  in  syngeneic  13762NF  rat  breast  carcinoma. 

Results;  A  significant  drop  was  found  in  the  mean  tumor  pOj  (23  to  9  mm  Hg,  p  <0.05)  within  90  min  after 
treatment  (30  mg/kg  of  combretastatin  A4  phosphate)  and  a  further  decrease  was  observed  at  2  h  (mean  2  mm 
Hg;  p  <0.01).  The  initial  changes  in  pOj  in  the  central  and  peripheral  regions  were  parallel,  but  by  24  h  after 
treatment,  a  significant  difference  was  apparent;  the  pOj  in  the  periphery  had  improved  significantly,  and  the 
center  remained  hypoxic.  These  data  are  consistent  with  DCE-MRI,  which  revealed  an  ~70%  decrease  in 
perfusion/permeability  (initial  area  under  signal-intensity  curve)  at  2  h  (p  <0.001).  The  initial  area  under 
signal-intensity  curve  recovered  fully  after  24  h  in  a  thin  peripheral  region,  but  not  in  the  tumor  center. 
Conclusion;  The  response  observed  by  DCE-MRI,  indicating  vascular  shutdown,  paralleled  the  pOj  measure¬ 
ments  as  expected,  but  quantitative  pOj  measurements  are  potentially  important  for  optimizing  the  therapeutic 
combination  of  vascular  targeting  agents  with  radiotherapy.  ©  2005  Elsevier  Inc. 

MRI,  Vascular  targeting  agent.  Vasculature,  Oxygenation,  Breast  tumor. 


INTRODUCTION 

Tumor  growth,  survival,  and  metastasis  depend  critically  on 
the  development  of  new  blood  vessels  (1).  Therefore,  ex¬ 
tensive  research  has  focused  on  developing  strategies  to 
attack  the  tumor  vasculature  (1,  2).  Tubulin-binding  agents 
(e.g.,  combretastatin  A4  phosphate  [CA4P]  and  ZD6126) 
represent  one  kind  of  vascular  targeting  agent  (VTA)  (3, 4). 
Promising  preclinical  studies  have  shown  that  such  agents 
selectively  cause  tumor  vascular  shutdown  and  subse¬ 
quently  trigger  a  cascade  of  tumor  cell  death  in  experimen¬ 
tal  tumors  (4, 5).  Although  massive  necrosis  can  be  induced, 
tumors  usually  regrow  from  a  thin  viable  rim.  Thus,  a 
combination  of  VTAs  with  additional  conventional  thera¬ 
peutic  approaches  will  be  required  (6,  7).  Several  studies 
involving  the  combination  of  VTAs  with  radiotherapy  (8- 
11)  or  chemotherapeutic  agents  (12)  have  shown  enhanced 
tumor  response. 

To  better  understand  the  mode  of  action,  and  hence, 
optimize  such  combinations,  in  vivo  imaging  approaches 
have  been  initiated  to  monitor  the  physiologic  changes 


resulting  from  VTA  administration  (13-15).  Dynamic  con¬ 
trast-enhanced  (DCE)-MRI  based  on  the  transport  proper¬ 
ties  of  gadolinium-diethylenetriamine  pentaacetic  acid  (Gd- 
DTPA)  is  the  most  commonly  used  imaging  approach  to 
study  tumor  vascular  perfusion  and  permeability.  DCE-MRI 
was  included  as  part  of  the  Phase  I  clinical  trials  of  CA4P 
(16,  17).  The  results  of  preclinical  and  clinical  DCE-MRI 
studies  have  shown  a  reversible  change  in  vascular  perfu¬ 
sion  in  the  tumor  periphery  after  a  single  dose  of  VTA 
(18-21).  For  combination  with  radiotherapy,  measurement 
of  tumor  oxygen  dynamics  will  be  especially  important, 
because  reduced  perfusion  can  induce  hypoxia,  potentially 
modulating  the  radiation  response.  A  number  of  studies 
have  reported  an  improved  response  when  administering 
VTAs  after  radiotherapy,  with  the  enhancement  reduced  or 
lost,  if  VTAs  were  administered  before  radiotherapy,  im¬ 
plying  increased  hypoxia  induced  by  VTAs  (8,  10).  Direct 
measurements  of  tissue  p02  using  the  Eppendorf  electrode, 
conducted  by  Horsman  et  al.  (22,  23),  found  increased 
hypoxia  3  h  after  CA4P  or  ZD6126  administration.  We  have 
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developed  a  method  for  measuring  tumor  oxygenation  and 
dynamics  based  on  nuclear  magnetic  resonance  echo 
planar  imaging  after  direct  intratumoral  injection  of  the 
reporter  molecule  hexafluorobenzene  (HFB)  called  “fluoro¬ 
carbon  relaxometry  using  echo  planar  imaging  for  dynamic 
oxygen  mapping”  (FREDOM)  (24,  25).  This  technique  pro¬ 
vides  PO2  measurements  at  multiple  specific  locations  si¬ 
multaneously  within  a  tumor  and  reveals  the  dynamic 
changes  at  individual  locations  with  respect  to  interven¬ 
tions.  We  have  previously  evaluated  the  tumor  oxygen 
response  to  varying  interventions  such  as  hyperoxic  gas 
breathing  (26,  27).  We  also  had  an  anecdotal  example  of 
PO2  response  to  a  tumor-selective  infarcting  agent  (28).  We 
have  now  applied  both  DCE-MRI  and  FREDOM  to  evalu¬ 
ate  tumor  perfusion/permeability  and  oxygen  dynamics  in 
response  to  CA4P  in  conjunction  with  confirmatory  histo¬ 
logic  examination. 

METHODS  AND  MATERIALS 

Tumor  model 

Rat  mammary  carcinoma  13762NF  was  implanted  syngenei- 
cally  in  a  skin  pedicle  surgically  created  on  the  fore  back  of  Fisher 
344  adult  female  rats  (n  =  25,  ~150  g,  Harlan),  as  previously 
described  in  detail  (29).  Of  the  25  rats,  9  were  used  for  DCE-MRI, 
10  for  FREDOM,  and  6  for  histologic  study.  The  Institutional 
Animal  Care  and  Use  Committee  approved  the  investigations. 

Drug  preparation  and  dosing 

CA4P  was  provided  by  OXiGENE  (Waltham,  MA).  CA4P  was 
dissolved  in  0.9%  saline  at  a  concentration  of  30  mg/mL  before 
each  experiment.  A  single  dose  of  30  mg/kg  CA4P  was  chosen  for 
this  study,  because  it  is  considered  a  clinically  relevant  dose  (18). 

MRI  experiments 

When  tumors  reached  ~1  cm  diameter  (~0.6  cm®),  MRI  was 
performed  using  a  4.7  T  horizontal  bore  magnet  with  a  Varian 
Unity  Inova  system.  Each  rat  was  given  intraperitoneal  ketamine 
hydrochloride  (120  juL;  100  mg/mL,  Aveco,  Fort  Dodge,  lA)  as  a 
relaxant  and  maintained  under  general  anesthesia  (air  and  1% 
isoflurane,  Baxter  International,  Deerfield,  IL).  A  27-gauge  butter¬ 
fly  catheter  (Abbott  Laboratories,  Abbott  Park,  IL)  was  placed 
intraperitoneally  for  infusion  of  CA4P  or  saline  alone.  For  DCE- 
MRI,  a  tail  vein  was  catheterized  using  a  second  27-gauge  butterfly 
catheter  for  contrast  agent  administration.  For  oximetry,  hexaflu¬ 
orobenzene  (50  /liL,  Lancaster,  Gainesville,  FL)  was  injected  di¬ 
rectly  into  the  tumor  along  two  or  three  tracks  in  a  single  central 
plane  of  the  tumor,  coronal  to  the  rat’s  body  using  a  Hamilton 
syringe  (Reno,  NV)  with  a  custom-made,  fine,  sharp  needle  (32- 
gauge),  as  previously  described  in  detail  (25).  A  tunable  (’H/'^T) 
volume  radiofrequency  (RF)  coil  was  placed  around  the  tumor¬ 
bearing  pedicle.  Each  animal  was  placed  on  its  side  in  the  magnet 
with  no  change  in  position  during  the  whole  study,  so  that  indi¬ 
vidual  regions  could  be  tracked.  A  thermal  blanket  was  used  to 
maintain  body  temperature. 

DCE-MRI 

Nine  tumor-bearing  rats  were  studied  before  CA4P  injection  (n 
=  6)  or  saline  alone  (n  =  3)  and  2  and  24  h  after  treatment.  On 


each  occasion,  a  series  of  T,-weighted  spin  echo  images  (TR  160 
ms,  TE  16  ms,  field  of  view  40  X  40  mm,  matrix  128  X  128,  voxel 
size  2.0  X  0.3  X  0.3  mm,  total  time  for  3  slices  23  s)  was  acquired 
before  and  after  a  bolus  injection  of  Gd-DTPA-BMA  (injection 
within  1  s;  0.1  mmol/kg.  Omniscan,  Amersham  Health,  Princeton, 
NJ)  on  three  2-mm-thick  cross-sections  parallel  to  the  animal.  Data 
were  processed  on  a  voxel-by-voxel  basis  using  software  written 
by  us  using  Interactive  Data  Language  (IDL),  version  5. 3/5 .4 
(Research  Systems,  Boulder,  CO).  For  each  slice,  the  tumor  was 
separated  into  central  and  peripheral  regions.  The  tumor  periphery 
was  taken  to  be  a  1 -2-mm-thick  rim  aligned  around  the  whole 
tumor.  Signal  intensity  vs.  time  curves  were  plotted  and  relative 
signal  intensity  changes  (AS/)  of  each  tumor  voxel  were  analyzed 
using  the  equation:  (AS/)  =  (S/^  —  S/*)/S/i,  where  S/^  refers  to  the 
enhanced  signal  intensity  in  the  voxel  and  S/*  is  defined  as  the 
average  of  the  baseline  images.  The  area  under  the  normalized 
signal  intensity-time  curve  (lAUC)  for  the  first  1.5  min  after 
Gd-DTPA-BMA  injection  was  integrated. 

tumor  oximetry — FREDOM 

A  separate  cohort  of  10  tumors  (7  treated  and  3  controls)  was 
used  for  pOj  measurement.  A  single  4-mm  slice  parallel  to  the  rat 
body  containing  the  strongest  fluorine  signal  was  chosen  for  the 
‘®F  MRI  PO2  studies.  *H  and  ‘®F  MR  images  were  acquired  using 
a  spin-echo  sequence.  Overlaying  the  ‘®F  MR  image  on  the  cor¬ 
responding  *H  image  revealed  the  distribution  of  HFB.  After 
conventional  MRI,  tumor  oxygenation  was  estimated  on  the  basis 
of  '®F  pulse  burst  saturation  recovery  echo  planar  imaging  relax¬ 
ometry  of  the  HFB,  as  previously  described  (24).  This  approach 
provided  pOj  maps  with  1.25  mm  in  plane  resolution  and  6-ju.L 
voxel  size  in  6.5  min.  The  spin-lattice  relaxation  rate  [R1  (s“‘)  = 
1/Tl]  was  estimated  on  a  voxel-by-voxel  basis  using  a  three- 
parameter  monoexponential  function.  p02  was  estimated  using  the 
relationship  pOz  mm  Hg  =  (R1  -  0.0835)/0.001876  (24).  The 
data  are  presented  in  bins  of  5  mm  Hg,  except  for  the  highest  and 
lowest  bins,  which  were  open  ended.  Before  CA4P  or  saline 
injection,  a  series  of  PO2  maps  was  acquired  with  respect  to 
respiratory  challenge  with  oxygen:  typically,  two  baseline  mea¬ 
surements,  three  with  oxygen  and  four  on  return  to  air.  Immedi¬ 
ately  after  the  last  (fourth)  air  measurement,  CA4P  (30  mg/kg)  or 
saline  (0.15  mL)  was  injected  intraperitoneally.  An  additional 
series  of  pOj  maps  was  acquired  after  10,  30,  60,  90,  and  120  min, 
and  finally  another  three  maps  while  breathing  oxygen.  The  24  h 
follow-up  study  comprised  two  measurements  with  air  and  four 
with  oxygen.  The  oxygen  challenge  was  included  to  evaluate 
vascular  function. 

Markers  of  vascular  perfusion  and  endothelium 

Six  animals  were  used  to  study  the  total  blood  vessels  and 
perfused  vessels  before  (n  =  2),  2  h  (n  =  2),  and  24  h  (n  =  2)  after 
CA4P.  The  blue  fluorescent  dye  Hoechst  33342  (Molecular 
Probes,  Eugene,  OR)  was  injected  into  the  tail  vein  of  the  anes¬ 
thetized  rats  at  a  coneentration  of  10  mg/kg  in  0.9%  saline  (0.1 
mL),  and  the  tumors  were  excised  1  min  later.  Tumor  specimens 
were  immediately  immersed  in  liquid  nitrogen  and  then  stored  at 
— 80'’C.  Immediately  after  cryostat  sectioning  (6  p.m  thick),  the 
slices  were  imaged  for  Hoechst  33342  under  ultraviolet  wave¬ 
length  (330-380  nm).  Perfused  vessels  were  determined  by  count¬ 
ing  the  total  number  of  structures  stained  by  Hoechst  33342  in  four 
fields  per  section  selected  to  show  high  perfusion  and  calculating 
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the  mean  number  of  vessels  per  millimeter  squared.  On  the  fol¬ 
lowing  day,  the  same  slices,  as  well  as  their  adjacent  slices,  were 
immunostained  for  the  endothelial  marker,  CD31.  A  primary 
mouse  anti-rat  CD31  monoclonal  antibody  (1:20  dilution,  Serotec, 
Raleigh,  NC)  was  added  and  incubated  for  2  h  at  37°C  in  a  humid 
box.  The  slides  were  incubated  with  Cy3-conjugated  goat  anti¬ 
mouse  secondary  antibody  (1:100  dilution,  Jackson  Immunore- 
search  Laboratories,  West  Grove,  PA)  for  1  h  at  37°C.  After 
mounting  with  Vectorshield  medium  (Vector  Laboratories,  Bur¬ 
lingame,  CA),  the  slides  were  observed  under  red  fluorescence 
(530-550  nm  excitation)  to  detect  anti-CD31  and  then  the  corre¬ 
sponding  Hoechst  33342  again  under  ultraviolet  light.  Image  anal¬ 
ysis  was  performed  using  Metaview  software  (Universal  Imaging, 
West  Chester,  PA). 


Statistical  analysis 

Statistical  significance  was  assessed  using  analysis  of  variance 
on  the  basis  of  Fisher’s  Protected  Least  Significant  Difference 
(StatView,  SAS  Institute,  Cary,  NC)  or  Student’s  t  tests. 


RESULTS 

DCE-MRI  findings 

Baseline  'H  MRI  showed  stable  signal  throughout  the 
tumor  (not  shown).  Immediately  after  a  bolus  injection  of 
Gd-DTPA-BMA,  the  signal  intensity  increased  signifi¬ 
cantly,  reaching  a  peak  after  23-69  s,  and  then  gradually 
decreasing  toward  baseline.  The  maximal  signal  enhance¬ 
ment  averaged  over  the  whole  tumor  in  all  three  slices 
ranged  from  ~52%  to  85%  for  the  nine  tumors. 

A  single  dose  of  CA4P  30  mg/kg  caused  a  dramatic 
decrease  in  signal  enhancement  at  the  2  h  point,  and  a  24  h 
follow-up  image  showed  signal  recovery  in  a  thin  tumor 
rim,  as  shown  for  a  representative  tumor  (Fig.  1).  Distribu¬ 
tion  of  the  lAUC  values  showed  a  distinct  shift  to  the  left  at 
both  2  and  24  h  after  treatment  (Fig.  lb).  The  highly 
enhancing  fraction  (lAUC  >1.5)  at  24  h  after  CA4P  reached 
the  pretreatment  level,  indicating  complete  recovery  in  the 
tumor  periphery  (Fig.  lb).  For  all  six  tumors  treated  with 


(b)  lAUC 


Fig.  1.  Dynamic  contrast-enhanced-magnetic  resonance  imaging  (DCE-MRI)  performed  before,  2  h,  and  24  h  after 
treatment  with  CA4P  (30  mg/kg).  Series  of  Tj-weighted  contrast-enhanced  images  acquired  at  these  times  with  respect 
to  infusion  of  bolus  injection  of  contrast  agent,  (a)  Normalized  images  of  signal  enhancement  23  s  after  contrast  agent 
administration  shown  for  representative  tumor.  Significantly  less  signal  enhancement  was  observed  for  whole  tumor 
region  2  h  after  treatment.  Full  recovery  is  apparent  in  tumor  rim  24  h  after  treatment,  (b)  Based  on  serial  images 
obtained  with  respect  to  contrast  agent  infusion.  Initial  area  under  signal-intensity  curve  (lAUC)  was  calculated  on 
voxel-by-voxel  basis.  Frequency  distribution  of  lAUC  values  is  shown.  Distributions  at  2  and  24  h  showed  significant 
increases  in  number  of  voxels  with  no  signal  enhancement  (lAUC  <0.05).  However,  frequency  of  highly  enhancing 
voxels  (lAUC  >1.5)  at  24  h  recovered  to  pretreatment  level. 
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CA4P,  a  significant  decrease  in  lAUC  was  observed  for  the 
whole  tumor  region  at  2  h:  the  tumor  peripheral  and  central 
region  was  reduced  by  65%  and  70%,  respectively  (p 
<0.05;  Table  1).  A  complete  recovery  in  lAUC  for  the 
tumor  periphery  was  seen  24  h  after  treatment,  but  the  tumor 
center  remained  low  (Table  1).  lAUC  frequency  data  from 
these  six  tumors  also  showed  a  significant  increase  in  the 
number  of  nonenhancing  voxels  (lAUC  <0.05),  and  the 
percentage  of  voxels  with  LAUC  >0.5  decreased  from  55% 
to  5%  2  h  after  administration  of  CA4P  (p  <0.05,  Table  2). 
In  control  tumors  treated  with  saline,  the  LAUC  was  similar 
at  all  points  (Tables  1  and  2). 

Tumor  oximetry — FREDOM 

Overlay  of  on  image  (Fig.  2a)  confirmed  that  HFB 
was  distributed  in  both  peripheral  and  central  regions  of  a 
central  plane  of  the  tumor.  In  the  series  of  echo  planar 
imaging  relaxation  data  sets,  typically  —40-100  voxels 
provided  an  R1  fit,  and  potential  pOj  value.  Because  noise 
may  give  an  apparent  relaxation  curve  (Rl)  fit,  data  were 
selected  by  applying  thresholds  of  T1  error  <2.5  s  and  T1 
error/Tl  ratio  <50%.  Only  those  voxels,  which  provided 
consistently  reliable  data  throughout  the  time  course,  were 
included  for  additional  analysis.  As  an  example  (Fig.  2b),  25 
voxels  qualified  for  all  these  prerequisites  and  could  be 
followed  through  the  17  measurements  from  baseline  to  2  h 
after  CA4P.  At  24  h,  31  voxels  were  traceable  through  the 
six  measurements  with  respiratory  challenge.  As  expected, 
on  the  basis  of  the  distribution  of  HFB  (Fig.  2a),  each  p02 
map  (Fig.  2b)  comprised  two  distinct  groups  of  voxels, 
representing  the  peripheral  and  central  location  of  HFB, 
respectively.  Variations  in  the  mean  pOj  with  respect  to 
intervention  for  these  two  groups  of  voxels  are  presented 
separately  in  Fig.  2c.  Before  CA4P  administration,  the  mean 
p02  was  45  mm  Hg  in  the  tumor  periphery  and  36  mm  Hg 
in  the  center  during  air  breathing.  Both  regions  responded 
significantly  to  oxygen  breathing  (p  <0.05),  and  this  was 
reversed  on  return  to  air  breathing.  A  significant  decrease  in 
PO2  was  detected  as  early  as  30  min  after  CA4P  adminis¬ 
tration  in  the  tumor  periphery  and  by  60  min  in  the  center  (p 


<0.05).  The  decline  continued  with  the  lowest  PO2  values 
(13  mm  Hg  and  10  mm  Hg)  at  the  2-h  point.  At  this  time, 
oxygen  challenge  no  longer  produced  an  increase  in  PO2  for 
either  region.  At  24  h  later,  the  peripheral  PO2  had  increased 
significantly  (24  mm  Hg),  but  was  still  significantly  lower 
than  the  pretreatment  value  (p  <0.05).  After  24  h,  PO2  again 
responded  significantly  to  breathing  oxygen  in  the  tumor 
periphery  (p  <0.01).  However,  PO2  in  the  central  regions  (4 
mm  Hg)  was  even  lower  at  24  h  than  at  2  h  and  did  not 
respond  to  oxygen  breathing.  Histograms  of  pooled  voxels 
assessed  for  the  seven  tumors  showed  significantly  de¬ 
creased  PO2  at  1,  2,  and  24  h  after  treatment  (p  <0.01).  The 
hypoxic  fraction  <5  mm  Hg  (HF5)  increased  significantly 
from  a  pretreatment  value  of  21%  to  36%  at  1  h  and  68%  at 
2  h,  but  was  only  34%  at  24  h  (Fig.  3).  For  the  seven  treated 
tumors,  two  tumors  showed  a  significant  decrease  in  p02  as 
early  as  30  min  after  treatment  (p  <0.05).  A  significant 
decrease  in  the  mean  p02  for  the  group  was  found  at  90  min 
after  CA4P  (23  ±  5  vs.  9  ±  3  mm  Hg;  p  <0.05)  and  a 
further  decrease  at  2  h  (2  ±  2  mm  Hg;  Table  3).  The  mean 
PO2  increased  significantly  to  15  ±  4  mm  Hg  at  24  h 
compared  with  2h  (p  <0.05).  For  the  control  tumors,  no 
significant  differences  in  p02  were  found  among  any  of  the 
measurements  (before  and  2  or  24  h  after). 

Histologic  examination  and  immunohistochemistry 

On  the  basis  of  the  perfusion  assessed  by  distribution  of 
Hoechst  33342,  a  significant  decrease  in  the  number  of 
perfused  vessels  was  evident  at  2  h  (21  ±  6/mm^  vs.  145  ± 
23/mm^  pre-CA4P,  p  <0.01),  followed  by  recovered  per¬ 
fusion  at  24  h  (114  ±  19/mm^).  Comparisons  of  perfused 
vessels  marked  by  Hoechst  33342  and  with  the  total  vessels 
labeled  by  anti-CD31  before  and  2  and  24  h  after  CA4P  are 
shown  in  Fig.  4. 

DISCUSSION 

In  common  with  previous  reports  in  the  literature  (15, 
18),  CA4P  caused  a  significant  reduction  in  tumor  perfusion 
within  2  h.  DCE-MRI  showed  that  the  lAUC  in  the  tumor 


Table  1.  Normalized  lAUC  by  DCE-MRI 


Group 

Mean  lAUC 

Baseline 

2  h 

24  h 

CA4P  (n  =  6)  (30  mg/kg) 

Periphery 

0.88  ±  0.06 

0.31  ±  0.08* 

0.86  ±  0.06 

Center 

0.44  ±  0.1  U 

0.13  ±  0.08* 

0.12  ±  0.05*** 

Control  (n  =  3)  (saline) 

Periphery 

1.09  ±0.04 

0.93  ±  0.13 

0.86  ±  0.05 

Center 

0.52  ±  0.19^ 

0.46  ±  0.05* 

0.44  ±  0.03* 

Abbreviations:  lAUC  =  initial  area  under  signal-intensity  curve;  DCE  =  dynamic  contrast 
enhanced;  CA4P  =  combretastatin  A4  phosphate. 

*  p  <0.05  from  baseline. 

^  p  <0.05  from  periphery. 

* p  <0.05  from  control. 
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Table  2.  Comparison  of  weakly  (lAUC  <0.05)  and  strongly 
(lAUC  >0.05)  responding  voxels  with  respect  to  DCE-MRI  and 


drug  administration 


Total  voxels  (%) 

Group 

lAUC 

Baseline 

2  h 

24  h 

CA4P  (n  =  6)  (30  mg/kg) 

<0.05 

16 

54*t 

41*t 

>0.5 

55 

5*t 

32*t 

Control  (n  =  3)  (saline) 

<0.05 

14 

15 

19 

>0.5 

62 

58 

51 

Individual  voxels  (total  =  10,508  in  CA4P  vs.  5376  in  control) 
categorized  according  to  lAUC. 

Abbreviations  as  in  Table  1. 

*  p  <0.01  from  baseline. 

^  p  <0.05  from  control. 

periphery  was  reduced  70%  2  h  after  a  single  dose  of  CA4P 
(30  mg/kg),  followed  by  a  full  recovery  by  24  h.  However, 
the  drop  in  the  lAUC  of  the  tumor  center  was  not  reversible 
(Table  1).  The  lAUC  reflects  tumor  blood  flow,  vascular 
permeability,  and  a  fraction  of  interstitial  space.  Histologic 
data  using  the  perfusion  marker,  Hoechst  33342,  confirmed 
a  significant  decrease  in  perfused  vessels  2  h  after  CA4P, 
with  recovery  evident  at  24  h  (Fig.  4). 

Vascular  shutdown  induced  by  VTA  is  expected  to  de¬ 
crease  tumor  PO2.  Horsman  et  al.  (23)  recently  reported  that 
tumor  PO2,  measured  by  Eppendorf  electrode,  was  signifi¬ 
cantly  decreased  3  h  after  ZD6126.  In  line  with  that  study, 
MRI  oximetry  showed  that  CA4P  caused  a  significant 
decrease  in  p02  within  90  min  and  a  further  reduction 
within  2  h.  As  for  perfusion,  p02  recovered  toward  pretreat¬ 
ment  levels  after  24  h  in  the  tumor  periphery  (Fig.  2C). 

Although  VTAs  cause  vascular  shutdown  and  lead  to 
extensive  hemorrhagic  tumor  necrosis,  rapid  tumor  re¬ 
growth  occurs  from  the  surviving  rim.  Many  preclinical 
studies  have  shown  that  a  single  dose  of  VTA  is  associated 
with  very  little  tumor  growth  delay  (4,  5,  23).  Thus,  a 
combination  with  other  therapeutic  modalities  (e.g.,  radio¬ 
therapy  and  conventional  chemotherapy)  will  be  needed  to 
achieve  complete  control.  The  rationale  for  such  combina¬ 
tions  is  that  VTAs  target  the  poorly  perfused  and  relatively 
hypoxic  regions  of  the  central  tumor,  and  conventional 
therapy  is  applied  to  attack  the  relatively  better  oxygenated, 
highly  proliferating  rim.  A  critical  issue  regarding  such 
combinations  is  the  timing  and  sequencing  of  the  agents. 
Thus,  a  primary  goal  of  this  study  was  to  monitor  longitu¬ 
dinal  PO2  changes  after  CA4P  to  facilitate  and  potentially 
optimize  such  combinations.  Others  have  demonstrated  en¬ 
hanced  therapeutic  efficacy  by  combining  CA4P,  ZD6126, 
or  5,6-dimethylxanthenone-4-acetic  acid  (DMXAA)  with 
radiotherapy  (8-11),  radioimmunotherapy  (30),  or  chemo¬ 
therapy  (12).  In  the  study  of  combined  ZD6126  with  radio¬ 
therapy,  Siemann  and  Rojiani  (8)  reported  increased  tumor 
cell  kill  by  administering  drug  immediately  (30  min  to  1  h) 
after  a  single  radiation  dose  compared  with  administering 
the  drug  24  h  or  1  h  before.  This  observation  coincides  with 
our  present  findings  that  tumor  p02  declined  1  h  after  drug 
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administration  and  had  not  recovered  fully  at  24  h  (Figs.  2C 
and  3). 

Oxygen  challenge  before  and  after  treatment  was  used  to 
compare  vascular  function  in  this  study.  In  common  with 
our  previous  observations  (31,  32),  oxygen  breathing  sig¬ 
nificantly  increased  tumor  oxygenation  in  mammary  carci¬ 
noma  13762NF  under  control  conditions  (Fig.  2  and  Table 
3).  The  PO2  response  to  CA4P  was  equally  effective  in  the 
peripheral  and  central  regions  through  2  h  after  administra¬ 
tion  (Fig.  2C).  PO2  decreased  in  both  regions  within  30  min 
and  continued  to  decrease  for  2  h.  Differential  behavior  with 
respect  to  oxygen  challenge  between  the  peripheral  and 
central  tumor  regions  was  observed  only  after  24  h.  At  24  h, 
PO2  in  the  tumor  center  was  even  lower  than  the  pretreat¬ 
ment  baseline  or  the  2-h  level  and  showed  no  response  to 
oxygen  inhalation  (Fig.  2).  However,  oxygen  breathing  pro¬ 
duced  a  significant  increase  in  PO2  of  the  peripheral  region, 
even  though  the  PO2  was  significantly  lower  than  before 
treatment  (Fig.  2).  These  observations  provide  further  evi¬ 
dence  that  CA4P  at  30  mg/kg  induces  irreversible  vascular 
damage  in  the  tumor  center,  and  peripheral  vessels  survived 
and  were  functional  by  24  h.  One  would  expect  that  hypoxia 
modifiers  (e.g.,  oxygen  breathing),  if  given  24  h  after  CA4P, 
might  improve  tumor  radiosensitivity. 

Hypoxia  in  solid  tumors  has  been  widely  recognized  as  a 
potent  factor  that  leads  to  resistance  to  radiotherapy  and 
some  anticancer  drugs  (33).  Recently,  increasing  evidence 
has  shown  that  tumor  malignant  progression  may  be  asso¬ 
ciated  with  a  hypoxic  microenvironment  (34,  35).  Given  the 
importance  of  oxygen,  many  techniques  for  monitoring  pOj 
have  been  developed  (25,  36).  Although  each  method  has 
specific  attributes,  many  are  highly  invasive  and  impractical 
for  longitudinal  studies  of  specific  regions  of  interest.  The 
Eppendorf  has  been  considered  the  gold  standard  for  p02 
measurement,  but  it  is  not  suitable  for  a  longitudinal  study 
of  the  type  reported  here.  In  vivo,  proton  MRI  (e.g.,  DCE- 
MRI  and  blood  oxygen  level-dependent  MRI)  provides  a 
noninvasive  approach  to  assess  tumor  vasculature,  particu¬ 
larly  useful  in  response  to  interventions,  but  neither  DCE 
nor  blood  oxygen  level-dependent  MRI  provides  a  straight¬ 
forward  correlation  with  p02  (37,  38).  FREDOM  not  only 
provides  PO2  values  simultaneously  at  multiple  specific 
locations  within  a  tumor,  but  also  reveals  dynamic  changes 
at  individual  locations  with  respect  to  interventions. 

Central  necrosis  develops  in  13762NF  tumors  even  at  a 
small  size  and  becomes  substantial  when  tumors  become 
bigger,  which  has  been  confirmed  previously  by  histologic 
examination  (32).  In  some  of  the  tumors  used  in  this  study, 
a  low-intensity  central  region  was  observed  in  Tj-weighted 
images,  in  which  little  enhancement  was  noted  after  infu¬ 
sion  of  the  contrast  agent.  After  treatment  with  CA4P,  the 
area  of  these  regions  of  low-signal  intensity  increased 
slightly  (data  not  shown).  A  similar  observation  was  re¬ 
ported  by  Beauregard  et  al.  (39)  in  human  colon  carcinoma 
xenograft.  In  common  with  our  previous  study  of  13762NF 
tumors  (31,  32),  we  found  considerable  intratumoral  heter¬ 
ogeneity  in  the  distribution  of  PO2  values,  ranging  from 
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Fig.  2.  (a)  Distribution  of  hexafluorobenzene  (HFB)  in  representative  tumor  (No.  4  in  Table  3).  Overlay  of  signal 
density  on  anatomic  image  indicates  HFB  in  both  peripheral  and  central  regions,  (b)  pOj  maps  obtained  from  same 
tumor  comprise  two  separated  groups  of  voxels  that  correspond  to  locations  of  HFB  on  anatomic  image;  25  individual 
voxels  were  traceable  from  pretreated  baseline  to  2  h  after  CA4P  with  oxygen  breathing;  and  31  voxels  could  be 
followed  24  h  after  treatment.  Significant  decrease  in  p02  was  evident  for  all  individual  voxels  after  CA4P.  p02  did  not 
respond  to  oxygen  inhalation  after  2  h;  24-h  maps  showed  improved  p02  and  significant  response  to  oxygen  breathing 
in  peripheral  region  (right),  but  not  in  central  region  (left),  (c)  Mean  p02  curves  shown  for  peripheral  (black  squares) 
and  central  (white  circles)  voxels  of  tumor.  *p  <0.05  from  baseline  air,  '^p  <0.05  from  24-h  air,  <0.05  from 
periphery. 


878 


I.  J.  Radiation  Oncology  •  Biology  •  Physics 


Volume  62,  Number  3,  2005 


0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 


Pre  (air) 

Mean  =  24±  2  torr 
Median  =  18  torr 


0.7 1 
0.6 


05 


0.4  i 


05 

05 


0.1 


CA4P  60  min 

Mean  =  14±  2  torr 
Median  =  13  torr 


Ik 


as  0.7 


CA4P  2  h 

Mean  =  3±  1  torr 
Median  =  2  torr 


15  25  35  45  55  65  75  65  95  >100 

pOj  (torr) 


0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 


14 


CA4P  24  h 

Mean  =  18±  2  torr 
Median  =  11  torr 


<5  15 


25  35  45  55  65  75  85  95  >100 

pOj  (torr) 


Fig.  3.  pOj  histograms  based  on  relative  numbers  of  voxels  in  seven  tumors  showing  left  shift  after  CA4P.  Ordinate 
ranges  (e.g.,  15  refers  to  10  mm  Hg  s  pOj  <  15  mm  Hg).  x  =  mean;  m  =  median. 


hypoxic  to  well  oxygenated.  Compared  with  our  previous 
data  on  13762NF  tumors,  we  found  the  average  baseline 
PO2  values  were  lower  in  the  present  study  (mean  25  mm 
Hg  vs.  63  mm  Hg).  This  difference  likely  resulted  from  an 
altered  anesthesia  protocol,  because  we  now  use  a  baseline 
of  air/isoflurane,  instead  of  33%  oxygen  in  nitrous  oxide 
with  metafane  (32).  However,  oxygen  challenge  produced  a 
remarkably  similar  extent  of  increase  in  p02  values  ob¬ 
served  in  both  studies  (72  mm  Hg  vs.  73  mm  Hg). 

Previous  studies  showed  little  macroscopic  redistribution 


of  HFB  during  a  period  of  hours  (40),  but  it  did  clear  from 
the  tumors  with  a  typical  half-life  about  600  min  (41).  The 
clearance  of  HFB  normally  precludes  chronic  long-term 
studies  of  oxygenation,  unless  additional  doses  of  HTO  are 
administered.  For  all  three  control  tumors  in  this  study, 
additional  HFB  (50  ju,L)  was  required  for  the  24-h  follow-up 
study.  However,  no  additional  HFB  was  required  for  studies 
24  h  after  CA4P,  presumably  owing  to  vascular  shutdown. 

These  results  provide  the  first  insight  into  regional  tumor 
oxygen  dynamics  in  response  to  CA4P  in  a  syngeneic  rat 


Table  3.  Tumor  oxygen  dynamics  assessed  by  ‘®F  MRI  with  respect  to  CA4P  treatment 


Group 

Case 

No. 

Size 

(cm^) 

PO2 

(mm  Hg) 

Baseline 

air 

Time  course  post  CA4P 

30  min 

60  min 

90  min 

2  h 

Oxygen 

24  h  base 

Oxygen 

CA4P  (30  mg/kg) 

1 

2.0 

36  ±  1 

NA 

NA 

NA 

2  ±  3* 

5  +  2* 

NA 

NA 

(n  =  7) 

2 

0.9 

12  ±  2 

NA 

NA 

NA 

-3  +  2* 

-3  ±  2* 

NA 

NA 

3 

1.2 

13  ±0 

35  ±  11 

11  ±  3 

9  ±  4 

1  ±  2* 

2  ±  2* 

5  ±  0* 

30  ±  6+ 

4 

0.7 

43  ±5 

29  ±  3* 

21  ±  3* 

20  ±  2* 

13  ±  1* 

16  ±  5* 

20  ±  2* 

48  ±  V 

5 

0.4 

12  +  2 

5  ±  3* 

9  +  3 

4  ±2* 

2  ±  1* 

7  ±  3* 

NA 

NA 

6 

0.7 

15  ±  1 

19  ±  3 

9  ±  2* 

4  +  2* 

1  ±  2* 

0  ±  3* 

12  ±  1 

73  ±  12+ 

7 

0.3 

28  ±  6 

29  ±  2 

16  ±4 

6  ±  2* 

1  ±  1* 

7  +  2* 

23  ±  2* 

79  +  6+ 

Mean 

0.9  ±  0.2 

23  ±  5 

23  ±  13 

13  +  5 

9  ±  3* 

2  +  2* 

5  ±  3* 

15  ±  4 

58  ±  11+ 

Saline  («  =  3) 

8 

0.6 

34  +  2 

47  +  5 

40  +  5 

36  ±  10 

31  ±  8 

88  ±  10* 

37  ±  0 

86  +  5+ 

9 

0.2 

47  ±  2 

62  ±  8 

65  ±  6 

64  ±  5 

61  ±  4 

114  ±  8* 

67  ±  1 

159  ±  7+ 

10 

0.6 

13  ±  1 

16  ±  6 

10  ±4 

9  ±  4 

12  ±  6 

35  ±  10* 

14  ±  1 

78  ±  15+ 

Mean 

0.5  ±  0.1 

31  ±  10 

42  +  14 

38  ±  16 

36  ±  16 

35  ±  14 

79  ±  23* 

39  ±  15 

108  ±  26+ 

Abbreviations:  NA  =  no  measurement;  other  abbreviations  as  in  Table  1. 
Data  presented  as  mean  ±  SE. 

*  p  <0.05  from  baseline. 

^  p  <0.05  from  24  h  air  (24  h  after  injection  of  CA4P). 
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Fig.  4.  Perfusion  marker  Hoechst  staining  before,  2  h,  and  24  h  after  CA4P  (left).  Vascular  endothelium  of  same  field 
immunostained  by  antl-CD31  (red,  middle).  Good  match  (right)  between  Hoechst  and  anti-CD31 -stained  vascular 
endothelium  found  in  pretreated  tumor.  Significant  reduction  in  perfused  vessels  detected  2  h  after  treatment  followed 
by  a  recovery  at  24  h.  Bar  =  100  /xm. 


tumor.  A  distinct  similarity  was  noted  between  the  results  of 
the  PO2  measurements  and  the  more  traditional  DCE,  but 
the  quantitative  PO2  values  provide  the  potential  for  exploit¬ 
ing  synergy  with  other  oxygen-dependent  therapies.  The 
observations  also  demonstrate  the  value  of  FREDOM  in 
assessing  dynamic  changes  in  regional  tumor  p02  in  vivo  in 


response  to  intervention.  We  believe  that  dynamic  measure¬ 
ments  are  particularly  valuable  for  understanding  the  mode 
of  action  of  therapeutic  response  to  VTAs.  Most  signifi¬ 
cantly,  these  measurements  lay  a  foundation  to  optimize  the 
timing  of  combination  therapy  involving  fractionated  radio¬ 
therapy  and  multiple  doses  of  VTAs. 
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ABSTRACT 

We  previously  suggested  that  the  two  time  constants  quantified  fi-om  the  increase  of  tumor  oxyhemoglobin 
concentration,  A[Hb02],  during  hyperoxic  gas  intervention  are  associated  with  two  blood  flow/perfusion  rates  in  well 
perfused  and  poorly  perfused  regions  of  tumors.  In  this  study,  our  hypothesis  is  that  when  cancer  therapy  is  applied  to  a 
tumor,  changes  in  blood  perfusion  will  occur  and  be  detected  by  the  NIRS.  For  experiments,  systemic  chemotherapy, 
cyclophosphamide  (CTX),  was  applied  to  two  groups  of  rats  bearing  syngeneic  13762NF  mammary  adenocarcinomas: 
one  group  received  a  single  high  dose  i.  p.  (200  mg/kg  CTX)  and  the  other  group  continuous  low  doses  (20  mg/kg  CTX 
i.  p.  for  10  days).  Time  courses  of  changes  in  tumor  A[Hb02]  were  measured  at  four  different  locations  on  the  breast 
tumors  non-invasively  with  an  Inhaled  gas  sequence  of  air-oxygen-air  before  and  after  CTX  administration.  Both  rat 
body  weight  and  tumor  volume  decreased  after  administration  of  high  dose  CTX,  but  continuous  low  doses  showed 
decrease  of  tumor  volume  only.  Baselines  (without  any  therapy)  intra-  and  inter-tumor  heterogeneity  of  vascular 
oxygenation  during  oxygen  inhalation  were  simitar  to  our  previous  observations.  After  CTX  treatment,  significant 
changes  in  vascular  hemodynamic  response  to  oxygen  inhalation  were  observed  fi-om  both  groups.  By  fitting  the 
increase  of  A[Hb02]  during  oxygen  inhalation,  we  have  obtained  changes  of  vascular  structure  ratio  and  also  of  perfusion 
rate  ratio  before  and  after  chemotherapy.  The  preliminary  results  suggest  that  cyclophosphamide  has  greatest  effect  on 
the  well  perfused  tumor  vasculature.  Overall,  our  study  supports  our  earlier  hypothesis,  proving  that  the  effects  of 
chemotherapy  in  tumor  may  be  monitored  non-invasively  by  using  NIRS  to  detect  changes  of  hemodynamics  induced 
with  respiratory  challenges. 

Keywords:  Breast  Cancer,  Cyclophosphamide,  Hemodynamics,  NIR  Spectroscopy,  Tumor  vascular  oxygenation 


1.  INTRODUCTION 

In  addition  to  surgical  resection,  many  other  types  of  cancer  therapy  are  available  for  patients  including  radiotherapy, 
photodynamic  therapy  and  chemotherapy.  Chemotherapy  plays  an  important  role  to  treat  cancers  even  though  it  has 
some  side  effects.  Currently,  the  effect  of  chemotherapy  is  monitored  by  MRJ  or  CT  that  can  measure  the  tumor  volume 
changes  during  cancer  treatment.  However,  it  can  take  up  to  3  weeks  to  detect  such  changes,  and  this  is  considerably  late 
for  clinicians  to  decide  whether  initial  therapeutic  strategy  should  be  continued  or  modified.  This  delay  in  detection  of 
chemotherapy  effect  can  reduce  the  quality  of  a  patient’s  life  and  ineffective  therapy  is  costly.  Therefore,  many 
researchers  are  trying  to  develop  tools  that  can  detect  the  early  response  to  cancer  treatment.  For  example,  Li  et  al.  have 
used  ^'P  nuclear  magnetic  resonance  spectroscopy  (NMRS)  to  measure  the  effectiveness  of  cyclophosphamide  (CTX) 
treatment  in  radiation -induced  fibrosarcoma  (RIF).*  They  found  that  the  ratio  of  inorganic  phosphate  to  other  phosphate 
metabolites  in  CTX  treated  group  was  significantly  decreased  during  the  tumor  growth  delay  period  compared  to  age- 
matched  controls.  Poptani  et  al.  studied  the  effects  of  CTX  treatment  in  RIF-1  tumors  in  terms  of  tumor  oxygenation 
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and  glycolytic  rate  changes  by  utilizing  MRS,  Eppendorf  electrode,  and  Redox  scanning.^  They  observed  that  CTX 
treatment  caused  reduction  in  glycolytic  rate,  a  significant  decrease  in  tumor  tissue  pOi,  and  also  an  increase  of  NADH 
levels  24  hours  after  the  treatment  while  tumor  volume  did  not  show  any  significant  difference  between  the  CTX-treated 
and  control  groups.  Zhao  et  al.  have  reported  significant  changes  in  rat  breast  tumor  perfusion  following  either  single 
dose  CTX  or  continuous  low  dose  “metronomic”  therapy.^ 

In  the  last  decade,  near  infrared  spectroscopy  (NIRS)  has  been  developed  to  examine  tissue  oxygenation  and  has  been 
widely  applied  to  investigate  hemoglobin  oxygenations  of  museles,'*’’’*  the  brain,’’*’’  and  animal  tumors.'®’' Since 
tumors  have  higher  vascular  density  and  also  higher  metabolism  than  normal  tissues,  total  herrioglobin  coneentration 
([Hbtotai]),  oxyhemoglobin  coneentration  ([Hb02]),  and  reduced  scattering  coefficient  (|is’)  were  used  as  markers  to 
identify  tumors  from  the  human  breast  by  using  NIRS.'"*’'*’'®  In  addition,  an  NIR  spectrometer  is  a  low-cost,  portable, 
and  real-time  display  instrument.  Therefore,  NIRS  has  a  good  potential  to  be  used  as  a  monitoring  tool  for  tumor 
treatment  planning  and  tumor  prognosis. 

We  have  previously  studied  breast  tumor  oxygenation  under  gas  intervention  using  NIRS  and  found  that  oxyhemoglobin 
concentration  changes  (A[Hb02])  during  gas  intervention  can  be  fitted  by  a  two-exponential  equation  containing  two 
time  constants.'’  Based  on  the  model,  we  formed  a  hypothesis  that  changes  in  oxygenated  hemoglobin  eoncentration 
result  from  well  perfused  and  poorly  perfused  regions  of  an  animal  tumor  to  explain  why  there  are  two  different  time 
constants  in  the  A[Hb02]  data.  The  model  further  allows  us  to  associate  the  signal  amplitudes  and  time  constants  to  the 
ratio  of  vascular  density  and  the  ratio  of  the  perfusion  rates  in  the  two  different  regions,  respectively. 

In  this  study,  we  applied  the  NIRS  system  to  monitor  the  tumor  oxygenation  changes  during  oxygen  intervention  before 
and  after  CTX  administration.  The  purpose  for  this  study  is  to  explore  the  NIRS  as  a  possible  tool  for  monitoring  tumor 
responses  to  chemotherapy.  This  work  is  based  on  the  following  hypothesis:  when  tumor  is  treated  with  chemotherapy, 
changes  in  blood  perfusion  and  vascular  density  in  the  tumor  will  occur  and  will  be  seen  as  changes  of  the  two  fitted 
parameters  from  the  NIRS  measurements.  In  addition,  by  developing  a  non-invasive  tool  for  monitoring  cancer  therapy, 
we  are  not  only  monitoring  the  reduction  of  tumor  size,  but  also  detecting  the  ehanges  of  tumor  physiological  conditions, 
which  are  essential  for  tumor  treatment  planning  and  tumor  prognosis. 


2.  MATERIALS  AND  METHODS 


2.1  Tumor  Model  and  Experimental  Procedure 

Rats  were  divided  into  three  groups  for  this  study.  Two  groups  were  treated  with  CTX  at  different  doses,  and  the  other 
group  was  administered  saline  instead  of  CTX  as  a  control  group.  Cyclophosphamide  was  chosen  as  a  chemotherapeutic 
agent  for  this  study  since  our  tumor  line  is  highly  responsive  to  all^lating  agents  and  platinum  chemotherapeutic 
agents.'’  CTX  is  an  antineoplastic  alkylating  agent,  and  it  has  been  used  to  treat  lymphomas,  cancers  of  the  ovary,  breast 
and  bladder,  and  chronic  lymphocytic  leukaemia.'*’"*  The  tumor  line  was  rat  mammary  adenocarcinomas  13762NF 
(cells  originally  provided  by  the  Division  of  Cancer  Therapeutics,  NCI),  and  the  tumors  were  implanted  in  the  hind  limb 
of  adult  female  Fisher  344  rats  (~200  g). 

The  rats  were  anesthetized  with  0.2  ml  ketamine  HCl  (100  mg/ml;  Aveco,  Fort  Dodge,  lA)  when  the  tumors  reached 
approximately  1  cm  in  diameter  and  maintained  under  general  gaseous  anesthesia  using  a  small  animal  anesthesia  unit 
with  air  (1  dm*/min)  and  1%  isoflurane  through  a  mask  placed  over  the  mouth  and  nose.  During  the  experiments,  the  rat 
was  placed  on  a  warm  blanket  to  maintain  body  temperature,  which  was  monitored  with  a  rectally  inserted  thermal  probe 
eonnected  to  a  digital  thermometer  (Digi-Sense,  model  91 100-50,  Cole-Parmer  Instrument  Company,  Vernon  Hills,  IL). 
Tumors  were  shaved  before  measurements  to  improve  optical  contact  for  transmitting  light.  A  pulse  oximeter  (model: 
8600V,  Nonin,  Inc.)  was  placed  on  the  hind  foot  to  monitor  arterial  oxygenation  (Sa02)  and  heart  rate.  For  the  single 
high  dose  group  (n=5),  a  light  source  and  four  detectors  from  a  multi-channel,  CW  (continuous  wave)  NIRS  (NIM,  Ine, 
Philadelphia,  PA)  were  attached  to  the  tumor  using  posts  and  swivel  post  clamps  (see  Figure  1(a)).  For  the  multi  low 
dose  (n=3)  and  control  (n=3)  groups,  we  have  used  four-channel,  frequency  domain  (FD),  NIRS  (ISS,  Champaign,  IL). 
In  the  latter  case,  the  four  sets  of  light  sources  replaced  four  detectors  shown  in  Figure  1(b),  and  one  detector  fiber  was 
placed  on  the  top  center  of  the  tumor  to  obtain  signals  from  four  different  regions  of  tumor. 
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Fig.  1.  Schematic  diagram  of  experimental  setup  for  the  multi¬ 
channel  NIR  spectroscopy  system,  (a)  A  temperature  probe  was 
located  rectally,  and  pulse  oximeter  was  placed  on  the  hind  paw  to 
monitor  the  arterial  oxygen  saturation,  (b)  A  light  source  probe  was 
placed  on  the  top  center  of  the  tumor,  and  four  detectors  were 
attached  on  four  sides  of  the  tumor  to  detect  the  regional  signals 
during  oxygen  intervention,  (multi-channel  CW  NIRS).  In  the 
setup  of  4  channel  FD  NIRS  system,  the  four  sets  of  light  sources 
replaced  four  detectors  and  one  detector  fiber  was  placed  on  the 
top  center  of  the  tumor  to  obtain  signals  from  four  different  regions 
of  tumor. 


All  the  measurements  were  performed  in  a  dark  room,  and  the  measurements  were  initiated  while  the  rats  breathed  air  for 
10  minutes  to  get  a  stable  baseline.  After  10  minutes  of  baseline  measurement,  the  inhaled  gas  was  switched  to  oxygen 
for  15  minutes,  and  then  back  to  air  for  15  minutes.  This  experimental  procedure  was  repeated  before  and  after  the  CTX 
treatment,  and  the  four  detectors  were  intended  to  be  located  at  the  same  positions  for  each  measurement  on  different 
days.  Using  an  ellipsoidal  approximation,  tumor  volume  V  (cm^)  was  estimated  as  V  =  (7t/6)  L  W  H  where  L,  W,  and  H 
are  the  three  respective  orthogonal  dimensions.  Raw  amplitude  data  from  four  locations  were  recorded  simultaneously 
during  the  experiments  and  processed  after  the  experiments  to  obtain  the  A[Hb02]  and  A[Hbtotai]-  The  amplitude  and 
time  constant  of  A[Hb02]  were  calculated  by  fitting  the  two-exponential  model  to  the  data  using  Kaleidagraph  (Synergy 
Software,  Reading,  PA).  The  corresponding  ratios  of  vasculature  coefficients  and  perfusion  rates,  i.e.,  Y1/Y2  fi/f2, 
were  also  calculated  to  show  static  and  dynamic  heterogeneities  of  the  tumors  at  different  locations. 

2.2  Measurement  System 

As  mentioned  above,  we  have  used  two  NIRS  systems  for  this  study;  a  multi-channel,  CW,  NIRS  system  was  used  in  the 
single  high  dose  study,  while  an  ISS,  4-channel,  frequency-domain,  NIRS  system  was  applied  in  the  continuous  low  dose 
study.  The  original  CW  NIRS  system  has  the  ability  to  measure  light  signals  from  eight  different  locations,  but  due  to  the 
finite  tumor  size,  we  used  only  four  detectors  to  monitor  tumor  vascular  oxygenation  dynamics  during  respiratory 
challenges.  In  common  with  our  previous  work’^’^®,  we  assume  that  oxyhemoglobin  and  deoxyhemoglobin  are  the  only 
significant  absorbing  materials  in  the  blood-perfiised  tumor  tissue  in  the  NIR  range.  The  absorption  coefficients 
comprise  the  extinction  coefficients  for  deoxyhemoglobin  and  oxyhemoglobin  multiplied  by  their  respective 
concentrations  (Eqs.  1  and  2)  at  730  nm  and  850  nm,  that  have  been  employed  in  our  multi-channel  NIRS  system. 

14"'“  =  +  eHb02’''’[Hb02],  (1) 

|ia*'‘’  =  eHb*'®[Hb]  +  eHb02''‘’[Hb02].  (2) 


5693-57  V.  4  (p.3  of  1 1)  /  Color:  No  /  Format:  Letter  /  Date:  2005-02-16  16:20:33 


SPIE  USE: 


DB  Check 


Prod  Check,  Notes: 


Please  verify  that  (1)  all  pages  are  present,  (2)  ail  figures  are  acceptable,  (3)  all  fonts  and  special  characters  are  correct,  and  (4)  all  text  and  figures  fit  within  the 
margin  lines  shown  on  this  review  document.  Return  to  your  MySPIE  ToDo  list  and  approve  or  disapprove  this  submission. 


Based  on  modified  Beer-Lambert’s  law,  the  data  presented  in  this  paper  were  analyzed  using  amplitude  values  to  find  the 
changes  in  absorption  (Eq.  3).  By  manipulating  Equations  1-3,  changes  in  oxygenated  hemoglobin,  deoxygenated 
hemoglobin  and  total  hemoglobin  concentrations  were  calculated  fi-om  the  transmitted  amplitude  of  the  light  through  the 
tumor  (Eqs.  4,  5  and  6). 


-  liaT  =  log  (Ag/Aj.)  /  L,  (3) 

7^0  8^0 

A[Hb02]  =  [-0.674*log  (Ag/A^)  +  1.1 17*log  "]  /  L,  (4) 

A[Hb]  =  [0.994*log  (Ag/A^)  -  0.376*log  (A^/A^)  ]/L,  (5) 

7'ift  ocn 

A[Hb,„, J  =  A[Hb]  +  A[HbO,]=[0.32*log  (Ag/A^)  +  0.741  *log  ”]  /  L,  (6) 


where  Ab  =  baseline  amplitude;  At  =  transition  amplitude;  L  =  optical  pathlengtb  between  source/detector.  The 
constants  contained  in  these  equations  were  computed  with  the  extinction  coefficients  for  oxy  and  deoxyhemoglobin  at 
the  two  wavelengths  used.^‘  Notice  that  these  coefficients  have  accounted  for  four  hemes  per  hemoglobin  molecule.  In 
principle,  L  should  be  equal  to  the  source-detector  separation,  d,  multiplied  by  a  differential  pathlength  factor  (DPF), 
i.e.,  L=d*DPF.  Little  is  known  about  DPF  for  tumors,  although  a  DPF  value  of  2.5  has  been  used  by  others.’'’  Since  our 
focus  is  on  dynamic  changes  and  relative  values  of  tumor  [Hb02]  in  response  to  oxygen  intervention,  we  have  taken  the 
approach  of  including  the  DPF  in  the  unit,  and  eq.  (4)  becomes  as  follows; 

730  R^rt 

A[Hb02]  =  [-0.674*log  (Ag/A,)  +  1.117*log  (A^/A^)  ]/d,  (7) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of  A[Hb02]  in  Eq.  (7)  is  mM/DPF. 

Since  the  wavelengths  of  light  sources  fi'om  the  ISS,  frequency-domain  system  were  750  nm  and  830  nm,  the 
corresponding  equations  for  A[Hb02],  A[Hb],  and  A[Hb(oiai]  are  modified  as  follows. 

750  830 

A[Hb02]  =  [-0.709*log  (Ag/A^)  +  1.404*log  (A^/A^)  ]/d,  (8) 

750  830 

A[Hb]  =  [0.9546*log  (AJ\)  -  0.4992*log  (A^/A^)  ]  /  d,  (9) 

A[Hb^^,J  =  A[Hb]  +  A[HbO,]=[0.2456*log  (Ag/A^f  *  +  0.9048*log  (Ag/A^)®'”]  /  d.  (10) 

2.3  Bi-exponential  Model  of  Tumor  Vascular  Oxygenation 

In  our  previous  report,’^  we  followed  an  approach  used  to  measure  regional  cerebral  blood  flow  (rCBF)  with  diffusible 
radiotracers,  as  originally  developed  by  Kety^^  in  the  1950’s.  By  applying  Tick’s  principle  and  defining  yas  the  ratio  of 
Hb02  concentration  changes  in  the  vascular  bed  to  that  in  veins,  we  arrived  at  Eq.  (11): 

AHbO2™““'^‘“‘=(0  =  ?H„[l-exp(-^/Y)]  =  A,[l-exp(-t/T)]  (1 1) 

where  yis  the  vasculature  coefficient  of  the  tumor,  H„  is  the  arterial  oxygenation  input  and /is  the  blood  perfusion  rate. 

If  a  tumor  has  two  distinct  perfusion  regions,  and  the  measured  signal  results  from  both  of  the  regions  (Figure  2),  then  it 
is  reasonable  to  include  two  different  blood  perfusion  rates,/  and/,  and  two  different  vasculature  coefficients,  y  and  y, 
in  the  model.  Therefore,  Eq.  (11)  can  be  modified  to  count  for  the  double  exponential  feature  observed  in  the 
experiments: 

AHb02™^™"“'’-(/)=  ^Ho[l-exp(-/t/;i)]+  ;iHo[l-exp(-///^)] 

=  Ai[l-exp(-//i:i)]  +  A2[l-exp(-f/T2)]  (12) 

where  /  and  y  are  the  blood  perfusion  rate  and  vasculature  coefficient  in  region  1  for  the  well  perfused  region, 
respectively;/  and  y  have  the  same  respective  meanings  in  region  2  for  the  poorly  perfused  region,  and  Ai  =  ;^Ho,  A2  = 
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^Ho,  Xi  =  jilfi,  t2  =  Then,  if  Ai,  A2,  Xi,  and  X2  are  determined  by  fitting  the  measurements  with  the  model,  we  can 
obtain  the  ratios  of  two  vasculature  coefficients  and  the  two  blood  perfusion  rates: 

=  /i  _  '^1/^2  ^13) 

Yi  A  ’  fi 

With  these  two  ratios,  we  are  able  to  understand  more  about  tumor  physiology,  such  as  tumor  vasculature  and  blood 
perfusion. 


3.  RESULTS 


3.1  Body  weight  and  tumor  volume  changes  during  chemotherapy 

Body  weight  and  tumor  volume  were  monitored  before  and  after  the  CTX  treatments  to  see  the  tumor  responses  and  side 
effects  fi-om  chemotherapy.  In  the  single  high  dose  treatment  group,  body  weight  decreased  until  6  days  after  the 
treatment,  but  later  increased  for  the  rest  of  days  of  observation.  Two  rats  among  five  in  this  group  failed  to  survive  at 
day  6  due  to  the  toxicity  from  the  high  dose  CTX  treatment.  Therefore,  the  data  shown  at  day  8  and  10  represent  the 
smaller  group  of  rats  which  survived  during  the  high  dose  treatment.  Tumor  volume  did  not  fiufher  decreased  after  day 
4.  (Fig.  2(b))  In  comparison,  rats  in  the  continuous  low  dose  group  initially  lost  weight  after  a  low  dose  of  CTX 
administration,  but  gradually  gained  the  weight  during  the  treatment,  presenting  low  toxicity  from  the  treatment.  This 
group  also  showed  a  significant  reduction  in  tumor  volume  during  the  treatment.  For  a  control  group,  saline  was  injected 
into  the  rats  instead  of  CTX,  and  a  gradual  decrease  of  body  weight  was  observed,  while  their  tumor  volumes  increased 
exponentially.  Changes  in  rat  body  weight  and  tumor  volume  were  normalized  to  day  0  (before  CTX  or  saline 
administration).  (Figures  2(a)  and  (b))  Solid  circles  represent  the  data  from  a  control  group,  and  open  squares  and  open 
diamonds  represent  the  continuous  low  dose  group  and  single  high  dose  group,  respectively. 


Fig.  2.  Normalized  changes  in  rat  body  weight  (a)  and  tumor  volume  (b)  during  the  saline  and  CTX  treatments. 


3.2  Intratumoral  heterogeneity  of  vascular  oxygenation  observed  by  the  multi-channel  MRS 

In  our  experiment,  we  have  utilized  either  one  light  source  and  four  detectors  and  from  the  CW  system  (NIM,  Inc)  or 
four  light  sources  and  one  detector  from  the  FD  system  (ISS,  Inc.)  for  the  rat  breast  tumor  measurements,  and  the  setups 
( Fig.  1).  After  10  minutes  of  baseline  measurement  with  air  breathing,  gas  was  switched  to  pure  oxygen,  causing  a  rapid 
increase  in  tumor  [Hb02].  These  changes  were  measured  simultaneously  from  four  locations  of  the  tumor.  Figure  3 
shows  a  representative  set  of  data  before  the  CTX  treatment,  with  the  DC  MRS  system.  Open  circles  show  the  raw  data 
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measured  by  multi-channel  NIRS,  and  the  solid  black  lines  represent  the  fitted  curves  using  our  bi-exponential  model  for 
hemodynamics  during  oxygen  intervention.'^  It  is  apparent  that  the  data  from  each  location  differs  though  there  are 
similar  trends  of  A[Hb02]. 
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Fig.  3.  Dynamic  changes  of  tumor  [Hb02]  from  four  locations  in  a  rat  breast  tumor.  The  rising  parts  of  A[Hb02]  from  the  four 
locations  were  fitted  using  a  double-exponential  expression.  Figures  3(a)-3(d)  were  taken  from  locations  #l-#4,  respectively.  In  this 
case,  the  tumor  was  not  treated  yet. 

Table  1.  Summary  of  vascular  oxygen  dynamics  determined  at  the  four  detectors  from  tumor  shown  in  the  Fig.  3. 


Location 

A, 

Az 

Tl 

■t2 

A1/A2  (—71/72) 

T1/X2 

Ml 

1  (Fig.  3a) 

0.013 

0.027 

0.48 

24.3 

0.48 

0.02 

24 

2  (Fig.  3b) 

0.029 

0.026 

0.09 

29.8 

1.12 

0.003 

373 

3  (Fig.  3c) 

0.008 

0.036 

0.15 

10.6 

0.22 

0.014 

16 

4  (Fig.  3d) 

0.021 

0.037 

0.12 

7.25 

0.57 

0.017 

34 

To  compare  the  data  taken  from  four  locations  of  the  tumor  more  clearly,  the  time  constants  and  amplitudes  from  the 
four  fitted  curves  are  summarized  in  Table  1.  The  ratios  of  Yi/72  and  fi/f2  characterize  tumor  vascular  structure  and  blood 
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perfusion  within  the  volume  of  tumor  interrogated  by  light.*^  In  principle,  when  Y1/Y2  is  close  to  1,  it  implies  that  the 
measured  optical  signal  results  equally  from  both  regions  1  (i.e.,  well  perfused  region)  and  2  (i.e.,  poorly  perfused 
region);  if  Y1/Y2  <  the  measured  signal  results  more  from  region  2  than  region  1  [Figures  3(a),  3(c)  and  3(d)].  As  Table 
1  demonstrates,  only  location  #2  has  a  ratio  of  Y1/Y2  slightly  higher  than  1,  and  the  readings  from  locations  #1,  #3  and  #4 
have  the  ratios  of  Y1/Y2  less  than  1.  This  may  suggest  that  the  tumor  volume  that  was  optically  interrogated  from  location 
#2  was  dominated  by  well  perfused  regions,  while  most  of  other  tumor  volumes  detected  from  locations  #1,  #3  and  #4 
are  composed  of  more  poorly  perfused  regions.  Furthermore,  all  the  ratios  of  /i/^  from  four  locations  of  the  tumor 
shown  in  Fig.  3  are  much  greater  than  1,  indicating  that  the  blood  perfiision  rate  in  well  perfused  region  is  much  greater 
than  that  in  poorly  perfused  region.  Especially,/]/^  from  location  #2  is  10  to  20  times  higher  than  those  from  locations 
#1,  #3,  and  #4,  showing  a  high  level  of  intratumoral  heterogeneity  in  dynamic  vascular  structure. 

3.3  Monitoring  vascuiar  hemodynamics  of  breast  tumors  before  and  after  chemotherapy 

The  tumor  hemodynamics  during  oxygen  intervention  were  measured  before  and  after  administration  of  CTX  and  saline. 
The  representative  data  from  the  control  group  and  continuous  low  dose  group  are  shown  in  Figs.  4(a)  and  4(b), 
respectively.  Similar  to  Figure  3,  open  symbols  are  the  raw  data  from  measurements,  and  solid  lines  are  the  fitted  curves 
using  our  double  exponential  model.  As  mentioned  before,  there  were  4  light  sources  placed  on  the  surface  of  tumor. 
Figure  4  shows  the  acute  and  then  gradual  changes  of  [Hb02]  after  switching  the  breathing  gas  from  air  to  oxygen,  and 
the  data  were  observed  at  the  same  (or  nearly  the  same)  location  of  tumor  from  day  0  to  day  6.  From  Figure  4(a),  we  can 
see  different  tumor  hemodynamics  at  different  days,  but  having  similar  trends  and  maximum  A[Hb02]  for  the  control 
group.  However,  the  data  taken  from  the  continuous  low  dose  group  show  quite  different  hemodynamics  throu^out  the 
treatment  days  (Fig.  4(b)).  Especially,  we  notice  that  the  fast  increase  part  became  much  smaller  at  Days  2  and  6 
compared  to  Day  0,  implying  a  significant  decrease  in  signal  from  the  well  perfused  region. 


Fig.  4.  Dynamic  changes  of  A??[Hb02]  taken  at  location  #1  from  a  rat  breast  tumor  before  and  after  administration  of  (a)  saline  and 
(b)  continuous  low  dose  of  CTX  (20  mg/kg  for  10  days).  The  rising  part  of  A[Hb02]  from  location  #1  was  fitted  using  the  double¬ 
exponential  expression. 

Figure  5  also  shows  the  changes  of  tumor  hemodynamics  during  oxygen  intervention,  before  and  after  a  single  high  dose 
of  CTX  treatment,  measured  by  the  CW  NIRS.  This  figure  clearly  demonstrates  that  we  can  observe  significant  changes 
in  tumor  hemodynamics  after  chemotherapy  by  using  respiratory  challenge  as  a  mediator.  The  fitted  parameters  from 
our  mathematical  model  are  summarized  in  Table  2  to  compare  the  changes  in  hemodynamic  parameters  before  and  after 
administration  of  CTX.  The  rising  part  of  A[Hb02]  from  location  #1  was  fitted  using  either  a  single-  or  double¬ 
exponential  expression. 
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Fig.  5.  Dynamic  changes  of  [Hb02]  taken  at  location  #1  from  a  rat  breast  tumor  before  and  after  a  single  high  dose  of  CTX  treatments 
(200mg/kg).  The  rising  part  of  A[Hb02]  from  location  #1  was  fitted  vising  either  a  single-  or  double-exponential  expression. 


Table  2.  Summary  of  vascular  oxygen  dynamics  determined  at  location  #1  from  the  tumor  shown  in  Fig.  5 

before  and  after  CTX  treatment. 


Day 

A, 

A2 

t2 

A1/A2  (=71/72) 

'C,/T2 

Ml 

0 

0.044 

0.031 

0.23 

25.21 

1.42 

0.0091 

156 

1 

0.032 

0.033 

0.13 

11.36 

0.97 

0.0114 

85 

3 

0.0087 

0.014 

0.089 

8.36 

0.62 

0.0106 

58 

5 

0.0087 

1.27 

4.  DISCUSSION 

Changes  of  rat  body  weight  and  tumor  volume  clearly  show  that  CTX  treatment  is  effective  for  the  tumor  type  that  we 
have  used  in  this  study.  For  the  control  group,  the  average  rat  body  weight  gradually  decreased  during  the  entire  course 
of  treatment,  which  implies  the  sickness  of  rats  possibly  due  to  the  tumor  growth  (cachexia).  (At  Day  10,  the  tumor 
volume  was  ~5  times  larger  than  that  at  Day  0.)  It  is  clear  that  there  is  a  different  effectiveness  of  CTX  treatment 
between  the  single  high  dose  group  (200mg/kg)  and  continuous  low  dose  group  (20mg/kg  for  10  days).  Bdth  of  the  CTX 
treatments  delayed  the  tumor  growth  and  even  further  reduced  the  tuinor  volume.  However,  a  single  high  dose  of  CTX 
treatment  caused  the  death  of  two  rats,  and  the  tumor  volume  was  not  decreased  further  4  days  after  the  treatment,  while 
the  continuous  low  dose  CTX  treatment  continued  to  provide  tumor  regression  without  causing  severe  sickness.  From 
this  observation,  it  is  obvious  that  continuous  low  dose  of  CTX  treatment  is  working  much  better  than  a  single  high  dose 
of  CTX  treatment  for  a  rat  mammary  adenocarcinomas  13762NF  tumor. 

NIRS  is  a  portable,  low  cost,  and  real  time  measurement  system  that  can  monitor  changes  of  vascular  oxygen  levels  by 
using  two  wavelengths.  We  have  previously  used  a  single-channel  NIRS  system  with  one  light  source  and  one  detector 
for  global  measurements  of  A[Hb02]  and  A[Hb,„,ai]  in  tumors  during  respiratory  challenges.’^’^®  Through  those 
experiments,  we  have  found  that  most  tumors  have  a  bi-phasic  behavior  in  A[Hb02]  increase  (i.e.,  a  rapid  increase 
followed  by  a  slow  and  gradual  increase)  after  switching  the  gas  fi-om  air  to  carbogen/oxygen.  To  explain  this  bi-phasic 
behavior,  we  developed  a  mathematical  model  and  formed  a  hypothesis  that  the  bi-phasic  behavior  of  A[Hb02]  during 
carbogen/oxygen  inhalation  results  from  two  different  vascular  regions  in  tumor  with  two  blood  perfusion  rates  and 
vascular  structures.'^ 
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By  giving  an  oxygen  intervention,  tumor  blood  vessels  are  acutely  subject  to  an  increase  of  [Hb02]  due  to  higher  supply 
of  oxygenated  blood  from  artery  compared  to  that  from  air  breathing.  However,  due  to  the  irregular  vascular  structure  in 
tumor,  the  well  perfused  regions  in  tumor  may  have  an  increase  in  [Hb02]  much  faster  than  other  parts  of  tumor  that  are 
poorly  perfused.  Therefore,  two  time  constants  obtained  from  tumor  hemodynamic  measurements  during  oxygen 
intervention  are  able  to  reveal  two  blood  flow/perfusion  rates  in  tumor,  more  precisely,  two  speeds  of  blood  flow  within 
the  tumor  blood  vessels.  More  recently,  we  have  shown  that  the  bi-phasic  increase  in  optical  density  changes  occurs 
when  there  exist  two  different  flow  rates  in  tumor  vascular  phantom.^^  The  bi-phasic  model  is  a  basis  of  our  current 
study  where  we  wish  to  detect  any  changes  in  vascular  structures,  hemodynamic  features,  or  perfusion  rate  within  a 
tumor  after  CTX  treatment. 

The  amplitude  and  time  constants  obtained  from  A[Hb02]  increase  (Fig.  5)  are  summarized  in  Table  2.  At  day  0,  we  can 
see  that  71/72  is  higher  than  1,  indicating  that  the  measured  signal  results  more  from  the  well  perftised  region  than  poorly 
perfused  region.  However,  this  ratio  becomes  less  than  1  after  injection  of  cyclophosphamide  (Day  1  and  3).  This  may 
be  explained  by  destruction  of  vascular  structure  in  tumor  after  chemotherapy.  We  expect  that  after  a  single  high  dose 
administration  of  CTX,  the  drug  circulates  in  the  blood  vessels  and  is  delivered  to  the  tumor  cells  more  in  the  well 
perfused  region  than  in  the  poorly  perfused  region.  This  will  lead  to  death  of  tumor  cells  in  the  well  perfused  region 
more  effectively  than  that  in  the  poorly  perfused  region,  eventually  resulting  in  decreases  in  tumor  volume  in  the  well 
perfused  region  more  than  in  the  poorly  perfused  region.  Then,  the  tumor  volume  containing  the  well  perfused  regions 
will  consequently  decrease,  so  will  the  contribution  of  detected  MR  signals  from  the  well  perfiised  region.  In  other 
words,  a  decrease  in  71/72  may  indicate  deceases  in  welt  perfiised  regions  in  tumor  volume,  after  the  administration  of 
CTX. 

As  shown  in  Table  2,  moreover,  the  perfusion  rate  ratio,  filfi,  was  also  decreased  after  a  single  high  dose  of  CTX 
administration.  At  Day  Q,f\lfi  was  very  high,  meaning  that  there  was  a  big  difference  of  perfusion  rate  between  the  well 
perfused  and  poorly  perfused  region  in  tumor.  However,  this  ratio  significantly  decreased  at  Day  1  and  3  after  CTX 
treatment,  representing  that  the  perfusion  rate  gap  between  the  well  perfused  region  and  poorly  perfiised  region  became 
much  smaller  than  that  at  day  0.  At  Day  5,  changes  in  [Hb02]  during  oxygen  intervention  do  not  show  any  bi-phasic 
behavior  anymore,  and  it  was  fitted  by  a  single-exponential  model.  This  may  indicate  that  most  of  tumor  cells  and/or 
tumor  vasculature  in  the  well  perfused  region  are  possibly  destroyed  by  the  effect  of  CTX,  resulting  in  quite  different 
hemodynamic  behavior. 


5.  CONCLUSION 

In  conclusion,  we  have  conducted  this  study  to  show  the  possibility  of  using  MRS  to  monitor  tumor  hemodynamics  in 
response  to  chemotherapy  by  comparing  the  changes  in  tumor  vascular  oxygenation  before  and  after  CTX  treatment. 
The  heterogeneity  of  tumor  vasculature  was  easily  observed  by  quantifying  the  blood  perfusion  rate  and  vascular 
coefficients  at  four  different  locations  of  the  tumor.  Tumor  hemodynamics  has  been  significantly  changed  before  and 
after  CTX  treatment  compared  to  the  saline-treated  control  group,  showing  high  possibility  of  the  MRS  system  to  be 
used  as  a  monitoring  tool  for  cancer  treatments.  Our  future  studies  will  include  the  development  of  MR  imaging  systems 
to  obtain  a  map  of  tumor  hemodynamic  changes  from  whole  tumor,  allowing  us  to  predict  the  efficacy  of  tumor 
treatment. 
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3799  Tumor  physiological  response  to  antivascular  agent  combretastatin  A4  phosphate 
assessed  by  magnetic  resonance  imaging 

■Dawen^  ^Zhaol,  Lan  Jiang,  Ammar  Adam,  Eric  W.  Hahn,  Ralph  P.  Mason.  UT  Southwestern 
Medical  Ctr.,  Dallas,  TX. 

Tumor  oxygenation  has  been  widely  recognized  as  a  potent  factor,  which  influences  tumor  response 
to  chemo-  and  radio-therapy.  The  vascular  targeting  agent,  combretastatin  A-4-phosphate  (CA4P) 
causes  tumor  vascular  shutdown  and  a  combination  of  CA4P  with  radiotherapy  has  shown  enhanced 
tumor  response.  We  have  recently  established  a  novel  magnetic  resonance  approach  to  measuring 
regional  tumor  oxygen  tension  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for 

Dynamic  Oxygen  Mapping)  with  hexafluorobenzene,  as  a  reporter  molecule.  Here,  the  effect  of 
CA4P  on  real-time  tumor  oxygenation  and  perfusion/permeability  was  monitored  in  vivo  in  NF 13762 

rat  breast  carcinomas  by  *^F  FREDOM and  ^H  dynamic  contrast  enhanced  (DCE)  magnetic 
resonance  imaging  (MRI).  Syngeneic  NF  13762  rat  breast  carcinomas  were  implanted  in  a  skin 
pedicle  on  the  foreback  of  Fisher  344  adult  female  rats.  When  tumors  reached  ~  1  cm  diameter  (~  0.6 
cm^)  MRI  was  performed  using  a  4.7  T  horizontal  bore  magnet  with  a  Varian  Unity  Inova  system.  A 
series  of  p02  maps  acquired  at  10,  30,  60,  90,  and  120  minutes  after  i.p.  CA4P  (30  mg/kg; 

OXiGENE,  Inc.  Waltham,  MA)  showed  a  significant  decrease  in  p02  of  some  tumors  as  early  as  30 
min  post  injection.  For  all  tumors  there  was  a  significant  drop  in  tumor  p02  within  90  min  after 
treatment  (mean  baseline  p02  -  23  torr  to  9  torr,  p  <  0.05)  and  a  further  decrease  was  observed  at  2  h 
(mean  =  2  torr;  p  <  0.01).  Intriguingly,  the  initial  changes  in  p02  in  central  and  peripheral  regions 
were  parallel,  but  by  24  h  post  treatment  significant  difference  was  apparent:  the  p02  in  periphery 

improved  significantly,  while  center  remained  hypoxic.  These  data  are  consistent  with  DCE  MRI, 
which  revealed  a  -70%  decrease  in  perfusion/permeability  (initial  area  under  signal-intensity  curve 
(lAUC))  at  2  h  after  30  mg/kg  CA4P  (p  <  0.001).  lAUC  fully  recovered  in  a  thin  peripheral  region, 
but  not  the  tumor  center  24  h  post-treatment.  Vascular  perfusion  marker  Hoechst  33342  staining 
confirmed  that  the  number  of  perfused  vessels  decreased  significantly  at  2  h,  and  then  recovered  at 
24  h.  We  conclude  that  magnetic  resonance  imaging  can  monitor  in  vivo  dynamic  changes  in  tumor 
vasculature  and  oxygenation  following  VTA  administration.  These  results  provide  the  first  insight 
into  regional  tumor  oxygen  dynamics  in  response  to  CA4P  in  a  syngeneic  rat  tumor.  While  dynamic 
contrast  MRI  indicates  vascular  shut  down,  the  critical  p02  measurements  are  potentially  more 

important  for  optimizing  therapeutic  combination  of  VTAs  with  conventional  therapy,  in  particular, 
irradiation.  Supported  by  POD  Breast  Cancer  IDEA  Award  (DAMP  1703103631.  in  conjunction 
with  NCI  ROl  CA79515/EB002762  and  the  Cancer  Imaging  Program,  a  P20  Pre-ICMIC  CA86354. 
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Hypoxic  cells  in  tumors  has  long  been  recognized  as  a  significant  factor  influencing 
tumor  response  to  chemo-  and  radio-therapy.  The  vascular  targeting  agent,  combretastatin 
A-4-phosphate  (CA4P)  causes  tumor  vascular  shutdown  and  a  combination  of  CA4P  with 
radiotherapy  has  shown  enhanced  tumor  response.  For  the  potential  combination  with 
radiotherapy,  measurement  of  tumor  oxygen  dynamics  will  be  especially  important.  Here, 
we  applied  a  novel  magnetic  resonance  imaging  (MRI)  approach  to  measuring  regional 
tumor  oxygen  tension:  FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging 
for  Dynamic  Oxygen  Mapping)  with  hexafluorobenzene,  as  a  reporter  molecule. 
FREDOM  allowed  us  to  quantitatively  determine  the  in  situ,  in  vivo  dynamic  action  of 
CA4P  on  real-time  tumor  oxygenation. 

Syngeneic  NF 13762  rat  breast  carcinomas  were  implanted  in  a  skin  pedicle  on  the 
foreback  of  Fisher  344  adult  female  rats.  When  tumors  reached  ~  1  cm  diameter  (~  0.6 
cm3)  MRI  was  performed  using  a  4.7  T  horizontal  bore  magnet  with  a  Varian  Unity 
Inova  system.  A  series  of  p02  maps  acquired  at  10,  30,  60,  90,  and  120  minutes  after  i.p. 
CA4P  (30  mg/kg;  OXiGENE,  Inc.  Waltham,  MA)  showed  a  significant  decrease  in  p02 
of  some  tumors  as  early  as  30  min  post  injection.  For  all  tumors  there  was  a  significant 
drop  in  tumor  p02  within  90  min  after  treatment  (mean  baseline  p02  =  23  torr  to  9  torr,  p 
<  0.05)  and  a  further  decrease  was  observed  at  2  h  (mean  =  2  torr;  p  <  0.0 1).  Intriguingly, 
the  initial  changes  in  p02  in  central  and  peripheral  regions  were  parallel,  but  by  24  h  post 
treatment  significant  difference  was  apparent:  the  p02  in  periphery  improved 
significantly,  while  center  remained  hypoxic.  In  contrast  to  tumor  behavior,  there  was  no 
significant  change  in  p02  of  thigh  muscles  of  rat  at  any  time  point  after  CA4P.  These 
data  are  consistent  with  dynamic  contrast  enhanced  (DCE)  MRI,  which  revealed  a  -70% 
decrease  in  perfusion/permeability  (initial  area  under  signal-intensity  curve  (lAUC))  at  2 
h  after  30  mg/kg  CA4P  (p  <  0.001).  lAUC  fully  recovered  in  a  thin  peripheral  region,  but 
not  the  tumor  center  24  h  post-treatment.  Vascular  perfusion  marker  Hoechst  33342 
staining  confirmed  that  the  number  of  perfused  vessels  decreased  significantly  at  2  h,  and 
then  recovered  at  24  h. 

These  results  provide  the  first  insight  into  regional  tumor  oxygen  dynamics  in  response  to 
CA4P  in  a  syngeneic  rat  tumor.  While  dynamic  contrast  enhanced  MRI  will  reveal 
vascular  changes,  dynamic  p02  measurements  are  potentially  more  important  for 
optimizing  combinations  of  VTAs  with  conventional  anti-cancer  therapies  against  breast 
cancers. 
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